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An extensive set of first-principles density functional theory calculations have been performed 
to study the behaviour of He, C and N solutes in austenite, dilute Fe-Cr-Ni austenitic alloys and 
Ni in order to investigate their influence on the microstructural evolution of austenitic steel alloys 
under irradiation. The results show that austenite behaves much like other fee metals and like Ni 
in particular. Strong similarities were also observed between austenite and ferrite. We find that 
interstitial He is most stable in the tetrahedral site and migrates with a low barrier energy of between 
0.1 and 0.2 eV. It binds strongly into clusters as well as overcoordinated lattice defects and forms 
highly stable He-vacancy (V m He„) clusters. Interstitial He clusters of sufficient size were shown to 
be unstable to self-interstitial emission and VHe n cluster formation. The binding of additional He 
and V to existing V m He n clusters increases with cluster size, leading to unbounded growth and He 
bubble formation. Clusters with n/m around 1.3 were found to be most stable with a dissociation 
energy of 2.8 eV for He and V release. Substitutional He migrates via the dissociative mechanism in a 
thermal vacancy population but can migrate via the vacancy mechanism in irradiated environments 
as a stable V2He complex. Both C and N are most stable octahedrally and exhibit migration energies 
in the range from 1.3 to 1.6 eV. Interactions between pairs of these solutes are either repulsive or 
negligible. A vacancy can stably bind up to two C or N atoms with binding energies per solute 
atom up to 0.4 eV for C and up to 0.6 eV for N. Calculations in Ni, however, show that this may 
not result in vacancy trapping as VC and VN complexes can migrate cooperatively with barrier 
energies comparable to the isolated vacancy. This should also lead to enhanced C and N mobility in 
irradiated materials and may result in solute segregation to defect sinks. Binding to larger vacancy 
clusters is most stable near their surface and increases with cluster size. A binding energy of 0.1 
eV was observed for both C and N to a [001] self-interstitial dumbbell and is likely to increase with 
cluster size. On this basis, we would expect that, once mobile, Cottrell atmospheres of C and N will 
develop around dislocations and grain boundaries in austenitic steel alloys. 

PACS numbers: 61.72.-y,61.82.Bg,71.15.Mb,75.50.Bb,75.50.Cc,75.50.Ee 



I. INTRODUCTION 

Steel, in its many forms, is the primary structural ma- 
terial in current fission and fusion systems and will be 
so for the foreseeable future. Carbon (C) and nitrogen 
(N) are both commonly found in steel, either as impor- 
tant minor alloying elements or as low concentration im- 
purities. In body centred cubic (bec) a-iron (a-Fe), it 
has been shown experimentally that C interacts strongly 
with vacancy point defects and more weakly with self- 
interstitial defects^ and can form so called Cottrell at- 
mospheres around dislocations^, influencing yield prop- 
erties and leading to strain ageing of the material. First 
principles (ab initio) calculations, as summarised in a 
recent review by Becquart and Domain^, support these 
findings and demonstrate that N exhibits similarly strong 
interactions. As such, both of these elements have a sig- 
nificant influence on microstructural evolution in bec Fe, 
even down to very low concentrations, and a detailed un- 
derstanding of their interactions and dynamics in steels 
is worthy of development, more generally. 

Helium (He) is produced in significant quantities in 
the high neutron-irradiation fluxes typically experienced 
by the internal components of fission reactors and in the 
structural materials for fusion systems by (n,a) transmu- 



tation reactions. In combination with the primary point 
defect damage typical of irradiated environments, the 
presence of He plays a critical role in the microstructural 
evolution of these materials. As a result of its low solu- 
bility in metals, He becomes trapped in regions of excess 
volume, such as dislocations, grain boundaries and, most 
strongly, in vacancies and vacancy clusters^—. As such, 
it aids the nucleation, stabilisation and growth of voids 
(He bubbles), resulting in swelling of the materia l 10 ' 13 ~ — . 
The formation of He bubbles has also been implicated in 
high-temperature embrittlement of material s 10 ' 17 : 18 . It 
is therefore of critical importance to gain a deep under- 
standing of the behaviour of He in these materials and 
the part it plays in the underlying mechanisms of mi- 
crostructural evolution. 

First principles electronic structure calculations offer 
the most accurate means to develop an atomic level un- 
derstanding of the dynamics and interactions of solutes 
and point defects in solids. As such, they play a cen- 
tral role in the development of a theoretical understand- 
ing of the microstructural evolution of irradiated materi- 
als, as part of a multi-scale modelling approach, such as 
that used in the FP6 project, PERFECT^ an d the FP7 
project, PERFORM602 .. 

The behaviour and interactions of He in a num- 
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ber of bcc and face-centred cubic (fee) metals have 
been studied using density functional theory (DFT) 
techniques^ ' 19 ' 21 " — . This database of He kinetics and 
interactions is essential for the interpretation of complex 
experimental results, such as those present in thermal He 
desorption spectra. A case in point is the work of Ortiz et 
al£^, who have developed a rate theory model based on 
DFT calculations of the kinetics and interactions of point 
defects, He and C in bcc Fo 25 ' 31 ' 32 . The model success- 
fully reproduces and interprets the existing experimental 
desorption results^. It is interesting to note that agree- 
ment with experiment was only possible once the effects 
of C were included, even though only 150 at. ppm of C 
was necessary; this again indicates the sensitivity of the 
microstructural evolution to C concentration. To date, 
however, there have been no ab initio studies of He in 
austenite, that is fee 7-Fe, or austenitic FeCrNi alloys. 
This, primarily, is a result of the difficulty in describing 
the paramagnetic state of these materials. 

Ab initio calculations have also been used to exten- 
sively study G^MM and N 4 -^ in bcc Fe. These cal- 
culations show excellent agreement with experimentally 
verifiable parameters, such as the migration energy bar- 
rier for C diffusion, where ab initio values of 0.86 e V 33 ' 36 , 
0.87 eV— and 0.90 eV— are in good agreement with the 
experimental value of 0.87 eV— For N, an equally good 
agreement is seen for the migration barrier, where a value 
of 0.76 eV was found by ab initio calculations^ 4 - and a 
value of 0.78 eV was found experimentally 39 . Calcula- 
tions in austenite are, however, limited primarily to so- 
lute dissolution, diffusion and their influence on the elec- 
tronic structure and local geometr y 33 ' 40 " — although cal- 
culations of vacancy-C binding have been performed 43 . 

In this work we present a detailed study of the en- 
ergetics, kinetics and interactions of He, C and N so- 
lutes in model austenite and austenitic systems using 
DFT. A full treatment of paramagnetic austenite and 
FeCrNi austenitic alloys would naturally take into ac- 
count the magnetic and composition dependence of the 
variables under study and while ab initio techniques 
are now becoming available to model the paramagnetic 
stated and calculations in concentrated alloys are cer- 
tainly achievable 4 ^, their complexity precludes a broad 
study of all the necessary variables relevant for radiation 
damage modelling. Previous studies have, instead, ei- 
ther taken ferromagnetic (fm) fee nickel (Ni) as a model 
austenitic syste m 19 ' 22 ' 47 or modelled austenite using a 
small set of stable, magnetically ordered states, as in our 
previous work 4 ^. The advantage is that a more detailed 
study is possible, but the level of approximation involved 
is certainly not ideal and careful use should be made of 
the results obtained. Here, we follow the same approach 
used in our previous worki£, performing our calculations 
in the two most stable ordered magnetic states of fee Fe. 
In addition, we present and compare the results of cor- 
responding calculations in fm Ni in order to make more 
general conclusions in Fe-Ni based austenitic alloys. 

In section [II] we present the details of our calculations. 



We then proceed to present and discuss our results for He, 
C and N solutes in defect free austenite and dilute Fe-Cr- 
Ni austenitic alloys in section IIIII and their interactions 
with point defects and small vacancy clusters in section 
IIVI before making our conclusions. 



II. COMPUTATIONAL DETAILS 

The calculations presented in this paper have been per- 
formed using the plane wave DFT code, VAS P 49 ' 50 , in 
the generalised gradient approximation with exchange 
and correlation described by the parametrisation of 
Perdew and Wang^i and spin interpolation of the cor- 
relation potential provided by the improved Vosko-Wilk- 
Nusair scheme^. Standard projector augmented wave 
potential o 53 ' 54 supplied with VASP were used for Fe, He, 
C, N, Ni and Cr with 8, 2, 4, 5, 10 and 6 valence electrons 
respectively. First order (N = 1) Mcthfcsscl and Pax- 
ton smearing^ of the Fermi surface was used throughout 
with the smearing width, a, set to 0.2 eV to ensure that 
the error in the extrapolated energy of the system was 
less than 1 meV per atom. A 2 3 k-point Monkhorst-Pack 
grid was used to sample the Brillouin zone and a plane 
wave cutoff of 450 eV. 

All calculations used supercells of 256 (±1, ±2,...) 
atoms, with supercell dimensions held fixed at their equi- 
librium values and ionic positions free to relax. For 
the relaxation of single configurations, structures were 
deemed relaxed once the forces on all atoms had fallen 
below 0.01 eV/A. For the nudged elastic band^ (NEB) 
calculations used to determine migration barriers an en- 
ergy tolerance of 1 meV or better was used to control 
convergence. Spin-polarised calculations have been per- 
formed throughout this work with local magnetic mo- 
ments on atoms initialised to impose the magnetic state 
ordering but free to relax during the calculation. The 
relaxed local magnetic moments were determined by in- 
tegrating the spin density within spheres centred on the 
atoms. Sphere radii of 1.302, 0.635, 0.863, 0.741, 1.286 
and 1.323 A were used for Fe, He, C, N, Ni and Cr re- 
spectively. 

We have performed our calculations in both the face 
centred tetragonal (fct) antiferromagnetic single layer 
(afml) and double layer (afmD) collinear magnetic ref- 
erence states for austenitic Fe (at T=0 K), which we 
will refer to as afmD Fe and afml Fe in what follows, 
respectively, using the same methodology as our previ- 
ous work 4 ^. Both of these structures consist of (ferro- 
magnetic (001) fee planes, which we will refer to as mag- 
netic planes in what follows, but with opposite magnetic 
moments on adjacent planes in the afml state and an 
up, up, down, down ordering of moments in adjacent mag- 
netic planes in the afmD state. The fee fm and fct fm 
states were found to be structurally unstable and spon- 
taneously transformed upon addition of a whole range 
of defects and solutes 4 ^. The fee ferromagnetic high- 
spin (fm-HS) state was, however, found to be stable to 
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isotropic effects and we have performed a select few cal- 
culations in this state for comparison with other work in 
the literature^. We have also performed a number of cal- 
culations in fee fm Ni, which we will refer to, simply, as 
Ni in what follows, where these results were not available 
in the literature. We take the lattice parameters for afml 
Fe as a = 3.423 A and c = 3.658 A, those for afmD Fe as 
a = 3.447 A and c = 3.750 A and take a = 3.631 A for 
fm-HS Fe. Calculations in Ni have been performed with 
an equilibrium lattice parameter of a = 3.522A. The cor- 
responding magnitudes for the local magnetic moments 
in bulk, equilibrium afml, afmD and fm-HS Fe were de- 
termined as 1.50, 1.99 and 2.57 /Xb, respectively^ and a 
local moment of 0.59 fiB was found in bulk, equilibrium 
Ni. Convergence tests indicated that local moments were 
determined to a few hundredths of /lib- 

We use elastic constants for our reference states, as 
determined previously^, or determined here using the 
same techniques. For fm-HS Fe, we find C\\ — 40 GPa, 
C12 = 240 GPa and C44 = — 10 GPa, which clearly shows 
instability to shear strains and tetragonal deformations, 
as C = C\\ — C12 = —200 GPa. It is, however, stable 
to isotropic deformations as the bulk modulus, B = 187 
GPa, is positive. For Ni, we find C n = 272 GPa, C 12 = 
158 GPa and C 44 = 124 GPa, which gives C = 114 GPa 
and B = 196 GPa, and shows that this material is stable 
to any strain deformations. 

We have determined the solution enthalpy for carbon 
in Fe and Ni using diamond as a reference state. The 
diamond structure was determined using the same set- 
tings as our other calculations but with sufficient k-point 
sampling to ensure absolute convergence of the energy. 
We found a lattice parameter of a = 3.573 A, in good 
agreement with the standard experimental value. 

We define the formation energy, E{, of a configuration 
containing nx atoms of each element, X, relative to a set 
of reference states for each element using 

E { = E-J2n^, (1) 
x 

where E is the calculated energy of the configuration and 
Ex is the reference state energy for element, X. We take 
the reference energies for Fe, Ni and Cr to be the energies 
per atom in the bulk metal, that is Fe in either the afml, 
afmD or fm-HS states, as appropriate, Ni in its fee fm 
ground state and Cr in its bec antiferromagnetic (afm) 
ground state. Details of the Fe and Ni reference states 
are given above, whereas for Cr an equilibrium lattice 
parameter of 2.848 A was found with a corresponding 
local moment of magnitude 0.87 /iB- For He, C and N 
the reference states were taken to be the non-magnetic 
free atom, as calculated in VASP. 

In a similar manner, we define the formation volume at 
zero pressure, Vf, of a configuration relative to the bulk 
metal by 

Vf = V(0) - n bu i k V bu ik, (2) 



where V(0) is the volume of the configuration at zero 
pressure, nbuik is the number of bulk (solvent) metal 
atoms in the configuration and Vbuik is the volume per 
atom in the defect-free bulk metal, which we found to be 

11.138, 10.712, 11.970 and 10.918 A 3 in afmD Fe, afml 
Fe, fm-HS Fe and Ni, respectively. For our calculations, 
V(0) was determined by extrapolation from our calcula- 
tions at the fixed equilibrium volume using the residual 
pressure on the supercell and the bulk modulus for the 
defect-free metal. 

We define the binding energy between a set of n 
species, {Ai}, where a species can be a defect, solute, 
clusters of defects and solutes etc., as 

n 

Eb(A u A n ) = ~ ^(^ A n) (3) 

»=i 

where E{(Ai) is the formation energy of a configuration 
containing the single species, Ai, and Ef (Ax, A n ) is the 
formation energy of a configuration containing all of the 
species. With this definition an attractive interaction will 
correspond to a positive binding energy. One intuitive 
consequence of this definition is that the binding energy 
of a species, B, to an already existing cluster (or complex) 
of species, {Ai, . . . , A n }, which we collectively call C, is 
given by the simple formula, 

E b (B,C) = E b (B,A 1 ,...,A n )-E h (A 1 ,...,A n ). (4) 

This result will be particularly useful when we consider 
the additional binding of a vacancy or solute to an al- 
ready existing vacancy-solute complex. 

We have quantified a number of uncertainties in the 
formation and binding energies presented in this work. 
Test calculations were performed to determine the com- 
bined convergence error from our choice of k-point sam- 
pling and plane-wave cutoff energy. For interstitial C 
and N solutes in a defect-free lattice, formation energies 
were converged to less than 0.05 eV and formation en- 
ergy differences, such as migration energies, to less than 
0.03 eV. For interstitial He the convergence errors were 
half of those for C and N. For configurations containing 
vacancies or self-interstitial defects, formation energies 
were converged to 0.03 or 0.07 eV, respectively, while 
binding energies were converged to 0.01 eV, except for 
the binding of He to a vacancy, where the error was 0.03 
eV. 

The zero-point energy (ZPE) contributions to the for- 
mation energy, which can be significant for light solute 
atoms, have not been calculated in this work. We per- 
formed calculations of the ZPE for He, C and N solutes 
in a number of test sites in afmD and afml Fe, keeping 
the much heavier Fe atoms fixed, which is equivalent to 
assuming they have infinite mass2£. The results showed 
that the ZPE contributions were consistently around 0.10 
eV in all cases, which we take as an estimate of the ZPE 
error on the formation energies of configurations contain- 
ing C, N and He. The variation with site was, however, 
surprisingly low at 0.01 eV, which we take as an estimate 
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of the ZPE error in formation energy differences, binding 
energies and in the solution enthalpy for C, given that 
ZPE contribution in graphite is very similar to in FeM. 

Performing calculations in a fixed supercell of volume, 
V, results in a residual pressure, P, for which an Eshelby- 
type elastic correction to the total and, therefore, forma- 
tion energyS£& of E COTT - = -P 2 V/2B, can be applied. 
As such, E corr - also serves to indicate the likely finite- 
volume error. For many of the configurations considered 
here these corrections are negligible compared to other 
sources of error. Where they are significant, however, 
their relevance is discussed at the appropriate points in 
the text. 



III. SOLUTES IN THE DEFECT FREE LATTICE 

A. Single solutes 

The formation energies for substitutionally and inter- 
stitially sited He, C and N solutes in the sites shown in 
Fig. [1] are given in Table U for Fe and Table HI for Ni. 
We found that the Eshelby corrections were negligible 
for substitutional He but could be as high as 0.02 eV in 
magnitude for interstitial He and 0.04 eV for interstitial 
C and N. The corresponding uncertainties in formation 
energy differences were around half of these values. We 
discuss the results for He first, followed by those for C 
and N. 



1. He solute 

We found that He exhibits a large, positive formation 
energy in all sites but is most stable substitutionally, 
which is consistent with existing DFT studies of He in 
other bcc and fee metals^ ' 19 i 21 ~ — . The standard expla- 
nation is that, as a closed shell noble gas element, bond- 
ing interactions should be primarily repulsive, leading to 
insolubility and a preference for sites with the largest 
free volum o 21 ' 26 . This result distinguishes He from other 
small solutes, such as C and N, which are more stable 
interstitially but also distinguishes it from substitutional 
alloying elements, such as Ni and Cr with formation en- 
ergy differences between substitutional and interstitial 
sites in Fe of 3.0 eV and abov o 48 ' 59 . 

In Fe, the influence of substitutional He on the local 
magnetic moments of atoms in its first nearest neighbour 
(Inn) shell was found to be similar to those for a va- 
cancy, being generally enhanced relative to the bulk mo- 
ment and by up to 0.38 /zb here. This is similar to He in 
bcc Fo 25 i 26 . Indeed, we found that if the He atom was re- 
moved from the relaxed substitutional configuration with 
no further relaxation, the local Inn Fe moments changed 
by less than 0.03 /ie- In contrast to the vacancy, however, 
where Inn Fe were displaced inwards by 0.09 and 0.02 A 
in afmD and afml Fe, respectively^ 8 -, the respective dis- 
placements around a substitutional He were, on average, 




FIG. 1: Substitutional (0) and interstitial octahedral (1), 
tetrahedral (2-3) and crowdion (4-6) positions (in black) in 
afmD Fe. The Fe atoms are shown in white with arrows to 
indicate the local moments. Magnetic planes are included to 
aid visualisation. The afmD Fe state, which has the lowest 
symmetry, is shown to uniquely identify all distinct positions. 
In afml Fe and Ni, positions 2 and 3 are equivalent by sym- 
metry, as are 5 and 6. In Ni, position 4 is also equivalent to 
5 and 6. 

outwards by 0.02 and 0.04 A. This contrast can also be 
seen in the formation volumes, which were found to be 
0.74 Vbuik and 0.96 Vbuik for a vacancy, compared to 1. 17 
Vbuik and 1.38 Vbuik for substitutional He in afmD and 
afml Fe, respectively. Results in Ni were very similar to 
Fe, with enhanced moments in the Inn shell around both 
a vacancy and substitutional He, a contraction of 0.04 A 
in the Inn shell around a vacancy and an expansion of 
0.02 A around substitutional He. The formation volume 
for substitutional He, at 1.02 Vbuik, was again found to 
be greater than that for the vacancy, at 0.66 Vbuik- 

The large formation energy difference, of around 2 eV, 
between substitutional He and the underlying vacancy in 
Fe and Ni [see Section lTV] . which must be due to chemical 
interactions, may seem at odds with the relatively inert 
behaviour of He mentioned above. However, similar re- 
sults in bcc Fe have been reproduced using simple pair 
potential o 60 ' 61 , which demonstrates that such a large en- 
ergy difference, once distributed over Inn and 2nn bonds, 
is commensurate with the relatively small forces observed 
on the neighbouring Fe atoms around substitutional He. 

In Fe, interstitial He was found to be most stable in 
the tetrahedral (tetra) site, the octahedral (octa) site be- 
ing the next most stable and lying 0.206 and 0.059 eV 
higher in energy in the afmD and afml states, respec- 
tively. There is no consistent ordering of the octa and 
tetra sites in ab initio studies of other fee metals, with the 
octa site being most stable in Ag2i, A122 and P d 21 ' 24 and 
the tetra site being most stable in Cu^i and N i 19 ' 21 ' 22 , as 
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He 



C 



N 



Config. 

sub (0) 
octa (1) 
tetra uu (2) 
tetra ud (3) 
[110] crow. (4) 
[Oil] crow, uu (5) 
[Oil] crow, ud (6) 



afmD Fe 


afml Fe 


afmD Fe 


afml Fe 


afmD Fe 


afml Fe 


E t 


Et 


Et 


Et 


Et 


E f 


(AE { ) 


(AEi) 


(AEf) 


(AE { ) 


(AE { ) 


(AEt) 


4.024 

(-) 


4.185 

(-) 


-6.981 

(-) 


-6.244 
(-) 


rlx (other) 


-5.153 
(-) 


4.669 


5.026 


-8.797 


-8.856 


-8.602 


-8.621 


(0.206) 


(0.059) 


(0.000) 


(0.000) 


(0.000) 


(0.000) 


4.529 


as 


-6.535 


as 


-6.917 


as 


(0.066) 


tetra ud 


(2.261) 


tetra ud 


(1.685) 


tetra ud 


4.464 


4.967 


-6.644 


-6.272 


-7.044 


-6.737 


(0.000) 


(0.000) 


(2.153) 


(2.585) 


(1.558) 


(1.884) 


r l v {o\ 


5.271 


-6.764 


-6.412 




-6.006 


rlx [o ) 


(0.303) 


(2.033) 


(2.445) 


rlx ) 


(2.614) 


4.827 


as 


-7.354 


as 


-7.000 


as 


(0.364) 


[011] crow, ud 


(1.443) 


[011] crow, ud 


(1.602) 


[011] crow, ud 


4.802 


5.188 


-7.487 


-6.744 


-7.218 


-6.328 


(0.338) 


(0.221) 


(1.310) 


(2.113) 


(1.384) 


(2.293) 



TABLE I: Formation energies, Ef, in eV, for substitutionally and interstitially sited solute atoms, as shown in Fig. [T] The 
formation energies in bold are for the most stable states. For He, which is most stable substitutionally, the most stable interstitial 
site is also highlighted. The formation energy differences, AE{ (in brackets), to the most stable interstitial configurations are 
also given, in eV. Where the configuration was found to be unstable the configuration to which it relaxed is given. The 
substitutional N configuration in the fct afmD state relaxed to one with an octa N at Inn to a vacancy. 



Config. 


He 

E f 
(AEt) 


C 
Et 
(AEt) 


N 
Et 
(AEf) 


sub (0) 


3.185 


-5.386 


-4.562 


(-) 


(-) 


(-) 


octa (1) 


4.589 


-8.422 


-7.520 








(0.129) 


(0.000) 


(0.000) 


tetra (2-3) 


4.460 


-6.764 


-6.497 




(0.000) 


(1.659) 


(1.023) 


(110) crow. (4-6) 


4.651 


-6.795 


-5.970 








(0.191) 


(1.628) 


(1.550) 


TABLE II: Formation energies, Et, in eV, for substitutionally 
and interstitially sited He, C and N atoms in Ni. The layout 
and data content is as in Table [I] 



our results for Ni confirm. In both Fe and Ni, however, 
He favours the tetra site, which gives a strong indication 
that the tetra site will also be the most stable interstitial 
site in concentrated Fe-Ni based austenitic alloys. 

The other interstitial sites considered here lie no more 
than 0.364 eV above the tetra site, suggesting many low 
energy migration paths for interstitial He, that is, in the 
absence of any lattice defects that can act as strong traps. 



The bilayer structure in afmD Fe breaks the symmetry 
of the octa site and a He atom placed there was found 
to spontaneously relax in the [001] direction (as defined 
in Fig. [IJ, to between layers of the same spin by 0.55 A. 
It is, perhaps, surprising that in both afml Fe and Ni, 
an octa-sited He was also found to be unstable to small 
displacements in many directions. We present the results 
of these calculations in Table IIIII 

It is particularly clear in the afml Fe data that lower 
energy configurations were found along all of our test di- 
rections, with He relaxing to between 0.23 and 0.62 A 
from the symmetrical position. The picture is less clear 
in afmD Fe, where He was generally found to relax to 
the lowest local energy minimum but other metastablc 
positions were found. In Ni, the drop in energy is far- 
less pronounced than in Fe but is still present, with He 
relaxing to stable positions 0.29 A from the centre along 
(100) directions and 0.54 A along (110) directions. These 
configurations are important, certainly as intermediate 
states for the migration of interstitial He, but also as po- 
tential transition states and already suggest a low migra- 
tion energy barrier. We study these possibilities in detail 
in Section lill B[ For completeness, we also tested for the 
presence of stable off-centre positions for tetra-sited He 
but relaxation always returned He to the symmetrical 
position. 

The displacements of Inn Fe atoms around 
interstitially-sited He were, unsurprisingly, found to 
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afmD Fe 


afml Fc 


Ni 


Config. 


El Ar 
(AE f ) 


Et 
(AEf) 


Ar 


El Ar 
(AE { ) 


octa 


sym. 


rlx 
octa [001] 


5.208 
(0.241) 


0.00 


4 - 617 0.00 
(0.157) 


octa 


[100] 


rlx 
octa [001] 


5.105 
(0.138) 


0.39 


4.607 Q29 
(0.147) 


octa 


[001] 


4 ' 812 0.30 
(0.348) 


5.026 
(0.059) 


0.50 


as 

octa [100] 


octa 


[001] 


4 - 669 0.58 
(0.206) 


as 

octa [001] 


as 

octa [100] 


octa 


[110] 


rlx 
octa [001] 


5.079 
(0.112) 


0.54 


4.589 Q54 
(0.129) 


octa 


[011] 


4799 0.58 
(0.335) 


5.035 
(0.068) 


0.23 


as 

octa [110] 


octa 


[011] 


rlx 
octa [001] 


as 

octa [011] 


as 

octa [110] 


octa 


[111] 


rlx 
tetra ud 


5.029 
(0.062) 


0.62 


rlx 
tetra 


octa 


[111] 


rlx 
octa [001] 


as 

octa [111] 


as 

octa [111] 



TABLE III: Formation energies, Ef , and formation energy dif- 
ferences, AE{ , in eV, to the most stable tetra site (in brackets) 
for octa-sited He atoms in Fe. He is either sited symmetrically 
(sym.) or has been displaced off-centre, in which case the di- 
rection of the displacement is used to label the configuration 
and the displacement length after relaxation, Ar is given, in 
A. The symmetrical position is as shown in Fig. [T] and direc- 
tions determined from that point with the coordinate system 
shown. When no stable local energy minimum was found the 
state to which the configuration relaxed is given. 



be greater than for the substitutional site. A tetra-sited 
He in afml Fe displaced its neighbours by 0.23 A. In 
afmD Fe, displacements of 0.22 and 0.32 A were found 
for tetra uu and tetra ud -sited He, respectively. The 
magnetic moments on the Inn Fe atoms were quenched 
relative to the bulk moments by 0.24 /ib in afml Fe 
and by 0.16 for the tetra uu site in afmD Fe but 
enhanced by 0.15 hb for the tetra ud site. We attribute 
this difference to the greater free volume into which Inn 
Fe atoms around a tetra ud site may be displaced. We 
found formation volumes of 0.82 Vbuik and 0.99 Vbuik for 
tetra-sited He in afml Fe and tetra-ud-sited He in afmD 
Fe, respectively. Once again, results in Ni were similar 
to Fe, with a 0.24 A displacement and moment quench 
of 0.09 /iB in Inn Ni atoms around a tetra-sited He and 
a formation volume of 0.78 Vbuik- 

In the most stable octa configuration in Fe, the local 
geometry is complicated by the displacement of He from 
the symmetrical position. For that reason, we define a 



local unit cell surrounding the octa site using the posi- 
tions of its six Inn metal atoms, which lie at the centres 
of the cell faces, and report on the lattice parameters of 
that cell. In both afml and afmD Fe, the local lattice 
parameter along [100] and [010] directions, a inn , is in- 
creased by 0.31 A relative to the bulk equilibrium lattice, 
with the local lattice parameter along the [001] direction, 
ci nn , exhibiting an increase of 0.26 A in the afmD state 
and 0.29 A in the afml state. The local moment on the 
Inn Fe atom that He is displaced towards is significantly 
quenched by 1.04 and 0.41 /j,b in the afmD and afml 
states, respectively. In contrast, the other Inn moments 
are moderately enhanced by between 0.03 and 0.17 /zb- 
In Ni, the most stable off-centre octa position is along 
(110) directions from the symmetrical position. The re- 
sulting local unit cell, which exhibits a very slight shear, 
has c inn 7^ Qinn, with di n n increased by 0.31 A relative 
to bulk and ci nn by 0.28 A. Local Inn Ni moments were 
found to be quenched by between 0.02 and 0.08 ^b- 

These findings suggest that the relative stability of 
tetra over octa He, which is opposite to the order sug- 
gested by free volume argument o 21 ' 26 , may be best as- 
cribed to the relative ease with which a tetra He may 
lower its purely repulsive interactions with neighbouring 
atoms by local dilatation. To further investigate this hy- 
pothesis in Fe we split the formation energy for unrelaxed 
and relaxed substitutional, octa and tetra He configura- 
tions into three terms, in a similar manner to the work 
of Fu et al£l. The first is the formation energy, Ef cL , of 
any underlying, atomically relaxed, defects, e.g. a single 
vacancy for substitutional He. The second is the mechan- 
ical energy, i?™ och - , required to deform the Fe matrix con- 
taining those relaxed defects to the exact positions found 
in the configuration under study. The third is the energy 
change from chemical interactions, _Ef hcm -, upon inser- 
tion of the solute into its final position with no further 
relaxation. We also define the insertion energy, E 1 / 13 ', as 
the sum of £J™ ech - and Ef hcm - i.e. the formation energy 
for insertion of a solute into any position in a relaxed 
Fe matrix containing any relevant defects. We take the 
insertion energy as a more appropriate measure of site 
preference than the (total) formation energy, E(. The 
results are given in Table IIVI 

The substitutional site is clearly the most favoured, 
even in the unrelaxed state and by at least 2.25 eV once 
relaxed. In the unrelaxed lattice, an octa He is signif- 
icantly more stable than a tetra He, as expected from 
purely repulsive interactions given the relative proximity 
of Inn Fe in the two sites. Under relaxation the chemi- 
cal bonding energy is significantly reduced and to a far 
greater degree in the tetra site. The positive mechanical 
deformation energy is also greater for tetra He but the 
net result is still to stabilise tetra over octa He. These 
results clearly show that the relative stability of He in 
tetra and octa sites can be understood as resulting from 
a balance between the energy required for local dilatation 
of the Fe matrix coupled with a purely repulsive Fe-He 
interaction, which we suggest could be easily modelled 
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Config. 




^-ichem. 




afmD Fe + He 


sub, unrelaxed 


0.136 


2.150 


2.286 


sub, relaxed 


0.167 


2.045 


2.212 


tetra, unrelaxed 


0.000 


6.778 


6.778 


tetra, relaxed 


1.330 


3.134 


4.464 


octa, unrelaxed 


0.000 


5.804 


5.804 


octa, relaxed 


0.755 


3.914 


4.669 


afml Fe + He 


sub, unrelaxed 


0.023 


2.662 


2.685 


sub, relaxed 


0.155 


2.073 


2.228 


tetra, unrelaxed 


0.000 


6.774 


6.774 


tetra, relaxed 


0.999 


3.968 


4.967 


octa, unrelaxed 


0.000 


6.081 


6.081 


octa, relaxed 


0.855 


4.171 


5.026 



TABLE IV: Mechanical deformation energy, £ f mcc \ and 
chemical bonding energy, _Ej hcm ', contributions to the total 
formation energy, Ef, and the insertion energy, E{ ns ' , for un- 
relaxed and relaxed substitutional, tetra and octa He solute 
configurations in afmD and afml Fe, in eV. The most stable 
octa configuration was used in both states and the tetra ud 
configuration was used for the afmD state. 

using a simple pair potential. 

In bcc Fe, the relative stability of tetra over octa He 
has been explained as resulting from strong hybridisa- 
tion of He p-states with Fe d-state a 21 ' 26 . However, we 
do not find the evidence for such strong hybridisation to 
be convincing. We suggest that a repulsive non-bonding 
mechanism also applies to bcc Fe and explains the differ- 
ence in a much simpler manner. The magnetic and po- 
larisation effects discussed by Seletskaia^ 6 - and Zu^ 1 - are 
a simple consequence of these non-bonding interactions 
and not He p-state, Fe d-state hybridisation. Formation 
energy calculations 21,26 show that octa-sited He is higher 
in energy both before and after relaxation, despite the re- 
laxation energy for octa He being greater than for tetra 
He. This results, primarily, from the very short Inn Fe- 
He separations in the octa site when compared to those 
for the tetra site and the relative strengths of the result- 
ing repulsive interactions. The fact that purely repulsive 
pair potentials for Fe-He interactions in bcc Fe are capa- 
ble of reproducing the relative stabilit y 60 ' 61 gives further 
support to our claim. 



2. C and N solutes 

The results for C and N solutes [in TablesHlandHI] show 
that both elements clearly favour the octa interstitial site 
in both Fe and Ni. Experimental observations show this 
to be the preferred site for C in an Fe-13wt%Ni-lwt%C 
austenitic alloy^ 2 -. One exception worth comment is that 
of substitutional C in afmD Fe, for which the insertion 



energy, _Ej ns -, which as discussed for He provides a more 
appropriate measure of site preference, is comparable to 
that for octa C. On further inspection we found that, 
due to the asymmetries in the afmD state, the initially 
on-lattice C atom relaxed to 0.77 A from the lattice site. 
While this displacement is certainly significant, the C 
atom remains closer to the substitutional site than to an 
octa position at Inn to the (vacated) lattice site and has 
been named to reflect this difference. Relaxation of the 
substitutional N configuration also resulted in displace- 
ment away from the lattice site but convergence was to 
a configuration with the N atom in an octa site at Inn 
to a vacancy. We performed calculations to test for the 
presence of any stable off-centre octa configurations for 
C and N but none were found, in contrast to the results 
for He. 

We discuss the influence of octa C and N solutes on the 
local lattice geometry in an identical manner to octa-sited 
He, that is, using a\ nn and ci nn . The results are pre- 
sented in Table [V] for both Fe and Ni, including results 
in fm-HS Fe, which was shown to be mechanically unsta- 
ble in our previous work 48 , but not to the isotropic strain 
exerted locally by an octa-sited solute. We include this 
extra state here to compare with the work of Jiang and 
Carter—. It is immediately clear that the geometrical 
influence of octa C is rather similar to octa N, although 
with slightly smaller dilatations for N. Local expansion 
is observed in all our reference states, although the ex- 
pansion of c in afmD and afml Fe is much less than for 
a. As a result, the local c/a ratio is significantly reduced 
relative to the bulk material, to 1.02 around a C solute 
in both afmD and afml Fe, which is in good agreement 
with the 3% tetragonal distortion found by Boukhvalov 
et al^, and to 1.01 and 1.02 around an N solute in afmD 
and afml Fe, respectively. 

The magnetic influence of octa C and N solutes is, 
again, very similar with significant quenching of the lo- 
cal moments on Inn solvent atoms seen in all reference 
states, as expected for magnetic atoms under compres- 
sion. In both afmD and afml Fe the effect is most pro- 
nounced in those neighbours lying within the same mag- 
netic plane as the solute, which also show the most signif- 
icant displacement, resulting in a quench of 0.72(0.66) /ib 
for C(N) in afmD Fe and of 1.25(1.37) in afml Fe. In 
fm-HS Fe, Inn moments are quenched by 0.48(0.57) //b 
around C(N) and in Ni a quench of 0.42 /^b was observed 
for both C and N. 

In addition to this local influence, we have investigated 
the dependence of the lattice parameters of our reference 
states on the fractional atomic compositions, x c and x N 
for C and N, respectively. For low concentrations, as 
studied here, the lattice parameters change linearly as 
a function of the fractional composition 6 ^. In this case, 
quantities such as Aa/(ax l c ), where Aa is the difference 
between the lattice parameter with and without solute 
atoms present, are dimensionless constants that com- 
pletely specify the linear expansion. Our calculations 
have been performed in supercells at the equilibrium lat- 
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afinD Fe afml Fe fm-HS Fe Ni 

C N C N_ 

0.174 0.145 0.183 0.170 



0.072 0.034 0.243 0.263 







In 




In 




0.321 


0.303 


0.305 


0.276 


A « 

AClnn 


a aqa 
0.080 


a riA 
0.048 


a 1 k /i 
0.154 


0.12/ 




1.01b 


1.013 


1.023 


1 AOO 

1.02.5 


Aa/(ax t x ) 


0.266 


0.265 


0.341 


0.327 


Ac/(cx f x ) 


-0.057 


-0.044 


0.026 


0.042 


Aa e s./(a c ff.x c x ) 


0.158 


0.162 


0.236 


0.232 


v t /v hum 


0.53 


0.54 


0.78 


0.76 


77S0L 


0.323 




0.263 





0.214 0.102 0.73 0.79 

-0.164 0.697 



TABLE V: Lattice parameter differences (Aai nn and Aci nn , in A) between those for the unit cell surrounding octa C and N 
solutes (ainn and ci nn ) and the bulk equilibrium lattice parameters and the local Ci an /ainn ratio. Linear expansion coefficients 
(Aa/(ax { x ) and Ac/(c:Ex)) for the dependence of the lattice parameters on the fractional atomic solute composition, £ x , for 
solute X. For afmD and afml Fe, the linear expansion coefficient for an effective lattice parameter, defined by a a g, = (a 2 c) 1//3 is 
also given. Fractional formation volumes, Vf/Vbuik for octa-sited C and N solutes are given. The solution energy, -E|°q , taken to 
dissolve graphite in each of the reference states is given, in eV. For comparison, our calculations in bcc fm Fe give 

£ f s G = 0.700 

eV. 



tice parameters, so we determine the linear expansion co- 
efficients by extrapolating to zero stress using the resid- 
ual stress that builds up on the supercell upon addition 
of a solute and a knowledge of the elastic constants [see 
Section |TT]. In afrnD and afml Fe we have also calculated 
the linear expansion coefficients for an effective lattice 
parameter, a c ff. = (a 2 c) 1 ^ 3 , as a means to compare more 
directly with experiment. 

The results [in Table |V] show that local expansion 
around the solutes leads to a net expansion of the cell, 
overall. The afrnD state of Fe does, however, exhibit a 
small contraction in c and the afml state shows very little 
expansion in c, when compared to that for a. Once again, 
this shows that the addition of C and N acts to reduce 
the c/a ratio, bringing the lattice back toward perfect fee. 
In austenite, experimental results by Cheng et alM- and 
presented by Gavriljuk et al^ show that Aa/(ax l c ) lies 
between 0.199 and 0.210, with Aa/(ax l N ) being slightly 
greater at between 0.218 and 0.224. Our results in afmD 
and afml Fe are in broad agreement with these values 
but do not differentiate between C and N. Results in fm- 
HS Fe are significantly different from experiment, which 
again shows the unsuitability of this state for modelling 
austenite. It is interesting to note that our results for 
Ni are consistent with those for austenite and do exhibit 
more expansion due to N than for C. Experimental results 
in austenitic FeCrNi alloys^ are comparable to those for 
pure austenite and the general agreement with our re- 
sults strengthens the case for using afmD and afml Fe or 
using Ni as model systems for austenite and austenitic 
alloys. 

We have determined the solution energy at fixed equi- 
librium volume, Ef°Q , taken to dissolve graphite into our 
four bulk states [TablefV]. We have done this by calculat- 
ing the solution energy relative to diamond and then ap- 
plying the commonly accepted experimental energy dif- 
ference of 20 meV/atom between the cohesive energies of 



diamond and graphite at T = K— . For comparison, 
we have calculated Ef°Q for C in bcc fm Fe and find a 
value of 0.70 eV, which is in good agreement with the 
experimental value of between 0.60 and 0.78 eV found 
by Shumilov et alM^. 

The solution energies in all three Fe states are signifi- 
cantly lower than for bcc fm Fe. This is consistent with 
the relatively higher solubility of C in austenite than in 
ferrite and with the well known experimental result that 
C stabilises austenite over ferrite, as seen in the phase 
diagram. The effect is most pronounced in the fm-HS 
state, where the reaction is exothermic, in good agree- 
ment with previous DFT calculations^. In combination 
with the results discussed above, this implies that at suf- 
ficient concentrations, C will act to stabilise the fm-HS 
state over the others, just as was found for Ni in fee Fe-Ni 
alloys^ 6 -. The same conclusions follow for N by a direct 
comparison of the formation energies for octa sited N [see 
Table H], for which we found values of -8.252 and -9.018 
eV in bcc fm and fee fm-HS Fe, respectively. 

Experimental results for the solution enthalpy of C 
in austenite^ yield a value of Ef°Q = 0.37el^ at the 
concentration studied here, which agrees to within 0.1 
eV with our calculations in afmD and afml Fe but not 
with those for the fm-HS state and again supports their 
suitability as reference states for paramagnetic austen- 
ite. Our calculations in the ferromagnetic state for Ni 
are in good agreement with previous DFT calculations 
of Siegel and Hamilton 68 . However, as they report, this 
value is higher than those found experimentally in high- 
temperature, paramagnetic Ni, which lie between 0.42 
and 0.49 eV. It is worth noting that their calculations in 
non-magnetic (nm) Ni, which they use to model the para- 
magnetic state, underestimate the experimental range at 
between 0.2 and 0.35 eV. We conclude that the calculated 
solution enthalpy for C in Ni is particularly sensitive to 
the underlying magnetic state. 
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B. Solute migration 

As a first step in the calculation of migration energies 
for He, C and N solutes we investigated whether a 32 
atom cell would be sufficient for this purpose. To do this 
we recalculated the formation energies for substitutional 
and interstitial He and C in afmD Fe using a 32 atom cell. 
We compare these with our 256 atom cell calculations 
[Table H] in Table ED 



Config. 


He 




C 










32 atom 


Error 


32 atom 


Error 




Et 


Et 


Et 


Et 




(AE f ) 


(AEt) 


(AEt) 


(AEt) 



Sub (0) 
octa (1) 
tetra uu (2) 
tetra ud (3) 
[110] crow. (4) 
[Oil] crow, uu (5) 
[Oil] crow, ud (6) 



4.039 

(-) 
4.730 
(0.151) 

4.607 
(0.028) 

4.579 

(0.000) 



0.015 

(-) 
0.061 
(-0.055) 

0.078 
(-0.038) 

0.115 

(0.000) 



rlx (3) rlx (3) 



4.897 
(0.318) 

4.866 
(0.287) 



0.070 
(-0.046) 

0.064 
(-0.051) 



-6.911 
(-) 
-8.798 

(0.000) 

-6.395 
(2.403) 

-6.544 
(2.255) 

-6.396 
(2.402) 

-7.209 
(1.589) 

-7.346 
(1.452) 



0.070 
(-) 
-0.001 
(0.000) 

0.140 
(0.142) 

0.100 
(0.102) 

0.368 
(0.369) 

0.145 
(0.146) 

0.141 
(0.142) 



TABLE VI: Comparison between calculations in 32 atom and 
256 atom cells in afmD Fe of the formation energies, Et, in 
eV, for substitutional and interstitial He and C solutes and 
formation energy differences, AEt, in eV, to the most stable 
interstitial configurations, highlighted in bold. The layout 
and data content of each column is as in Table HI The column 
headed, 32 atom, contains the results for the 32 atom cell and 
the column headed, Error, contains the difference between the 
32 atom and 256 atom results, which we take as an estimate 
of the finite volume error in the 32 atom cell. 



effects. 

For C, the choice of cell size actually changes the rel- 
ative stability of the [110] crowdion and tetra ud config- 
urations. This is important as these two are transition 
states on two distinct migration paths for C (as will be 
shown in what follows). The small cell would, there- 
fore, give the wrong minimum energy path (MEP), as 
found previously for C in fm-HS Fe2^. Closer inspec- 
tion of the [110] crowdion configuration showed that the 
periodic boundary conditions in the smaller cell applied 
unphysical constraints on the displacements of Fe atoms 
at Inn to C and along the crowdion axis generally, which 
resulted in a significant buckling, moving the C atom to- 
wards the tetra uu site, that is along [001], by 0.71 A. In 
the larger cell these constraints are not present, result- 
ing in a significantly lower formation energy and while 
there is still a small displacement towards the tetra uu 
site of 0.18 A this is to be expected given the asymmetry 
present in the afmD state. 

As a final test, we investigated the case of C in fm-HS 
Fe, where Jiang and Carter have determined a migration 
barrier in a 32 atom cell^. They found that the (110) 
crowdion site is an intermediate site for C migration, ly- 
ing only 0.01 eV below the transition state energy and 
0.98 eV above the stable octa site. Our calculations in 
a 32 atom cell agree well with this finding, with an en- 
ergy difference of 1.01 eV between the (110) crowdion 
and octa sites for C. However, when we repeated the 
calculations in a 256 atom cell, we found that a configu- 
ration with C in the (110) crowdion was structurally un- 
stable and spontaneously transformed as a result of the 
non-isotropic stress on the Fe lattice. By contrast, the 
isotropic stress from an octa-sited C only led to local re- 
laxation of the Fe matrix and maintenance of the crystal 
structure. We conclude that the 32 atom cell effectively 
imposed artificial constraints that allowed a, seemingly, 
sensible migration barrier to be determined. 

Overall, we find that the finite size effects in the 32 
atom cell are too significant and while some intermedi- 
ate cell size between the two investigated here may be 
sufficient, we have performed our migration energy cal- 
culations in the 256 atom cell. 



There is a significant size effect on the formation ener- 
gies in the 32 atom cell, except for substitutional He and 
octa C, where the formation energies are within errors 
of those in the 256 atom cell. The formation energies 
are greater in the 32 atom cell, as expected from volume- 
elastic effects, by between 0.06 and 0.12 eV for interstitial 
He and by between 0.00 and 0.37 eV for interstitial C. 
Formation energy differences to the most stable inter- 
stitial configuration also exhibit a significant size effect, 
with the smaller cell underestimating them by between 
0.04 and 0.06 eV for He configurations and overestimat- 
ing them by between 0.10 and 0.37 eV for C. It is reason- 
able to assume that the migration energy, which is itself 
a formation energy difference, will suffer from similar size 



1. Interstitial He migration 

The migration of interstitial He is relevant in the initial 
stages after He production by transmutation, a-particle 
irradiation and at sufficiently high temperatures for He 
to escape from defect traps. The migration of He between 
adjacent tetra sites (that is, between sites at Inn on the 
cubic sub-lattice of tetra sites) can proceed along many 
distinct paths, with their corresponding transition states 
defining the energy barrier for the transition. A direct 
path would lead to an intermediate state with He in the 
crowdion position but the energy differences to the tetra 
configurations (in Tablesllland lllll) suggest that the direct 
path is not the MEP and that the transition state has He 
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FIG. 2: Possible migration paths for interstitial He in the 
afmD Fe lattice. Paths are shown for Inn jumps from ini- 
tial to final tetra positions (black circles) via off-centre octa 
intermediate transition state positions (grey circles). The Fe 
atoms (white circles) are shown with arrows to indicate the 
local moments. The symmetry of the afmD state leads two 
distinct tetra sites (uu and ud) and three distinct Inn jumps, 
as shown. In the afml state paths 1 and 3 are equivalent but 
still distinct from path 2. Coordinate axes are as in Fig. [1] 



0.4 



0.3 



> 



0.2 



3 



0.1 



1 


1 


1 


1 1 

\ o afmD Fe - 
<S □ afml Fe 


/ Path 1 \ 


/ Path 2 T CD Path 3 ffl 




i F^ m \ \ 1 L. 


T 




1 


i 



Tetra 
ud 



Tetra 



Tetra 

ud uu/ud 

Reaction Coordinate 



Tetra 
uu 



FIG. 3: Formation energy difference, AE{, to the lowest en- 
ergy tetra configuration along the distinct migration paths for 
interstitial He in Fe, as shown in Fig. [2] Positions of the tetra 
configurations and the intermediate configurations, L;, along 
path i are labelled. In the afml state, the data for path 3 has 
been omitted as it is equivalent to path 1. The arrows in- 
dicate the expected lowering of the migration barrier heights 
if a reorganisation of magnetic moments is allowed along the 
migration path. 



in an off-centre position. We show representative paths 
for the three distinct Inn jumps in afmD Fe in Fig. [5] 

We have performed NEB calculations for He migra- 
tion in Fe along these paths and show the formation 
energy difference to the most stable interstitial config- 
uration against a, suitably chosen, reaction coordinate in 



Path, i 


afmD Fe 


afml Fc 


Ni 


1 


0.335 


0.070 


0.129 


2 (ud to uu) 


0.349 


0.119 




2 (uu to ud) 


0.283 






3 


0.160 







TABLE VII: Migration energy barrier height, E m , in eV, for 
the migration of interstitial He along the distinct paths iden- 
tified in Fig. [2] In afmD Fe, path 2 is asymmetrical and 
the direction of migration has been identified by the initial 
and final tetra sites. In Ni, all paths are equivalent and the 
migration energy is given by the formation energy difference 
between the octa [110] and tetra He configurations from Ta- 
bles M and UHl 



Fig. [3l with the corresponding migration barrier heights 
given in Table IVIII It is immediately obvious that the 
results for the two Fe reference states differ significantly. 
This is not surprising, however, given that typical mag- 
netic effects can be of the same order of magnitude as the 
migration barrier height [see Table ITTT] . The high barriers 
along paths 1 and 2 in afmD Fe are because the lowest- 
energy tetra site is between layers of unlike moment and 
so not adjacent to the lowest energy octa intermediate 
site, which lies off-centre between like-spin layers [Ta- 
ble IIII] . A wholesale reorganisation of spins would lower 
the barriers along these paths and would be preferred in 
the paramagnetic state. This problem is not present for 
path 3, resulting in a significantly lower barrier, which is 
more consistent with those found in the afml state, where 
a more uniform distribution of energies around the octa 
site exists [see Table IIII] . Path 1 in the afml state and 
path 3 in the afmD state both show a double peaked 
structure with weakly stable octa [001] and octa [001] in- 
termediate states, respectively. These intermediates are 
equidistant from four tetra sites, resulting in the same 
energy barrier for Inn and 2nn jumps on the tetra sub- 
lattice. The same cannot be said for migration along path 
2, which proceeds via a (near-)octa [110] transition state 
in both reference states. In the afml state, there appears 
to be a very shallow minimum at I2, that is the off-centre 
octa [110] configuration, but with a depth of 0.007 eV this 
may well be just an artefact of the convergence criteria 
as it is less than the expected error for formation energy 
differences. The data also exhibits a shoulder between 
the tetra ud and I2 sites, which we suggest results from 
close proximity to the octa [111] configuration. It seems 
reasonable to suggest that the barriers for 2nn jumps that 
cross a magnetic plane will be close in energy to those for 
path 2, given the additional data in Table Mil Overall, 
our findings suggest an energy barrier for interstitial He 
migration in austenite that is below 0.35 eV and more 
likely in the region between 0.1 and 0.2 eV. 

Such low migration barriers are typical of all metals 
for which data is available. Ab initio calculations find a 
value of 0.10 eV for fee 0.07 eV for fee Pd^i and 

0.06 eV for bec Fe2^, W— and V—. Experimental vali- 
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elation of these results is not forthcoming, primarily due 
to the low temperatures involved and the complications 
of He interactions within the material. In bec W, Wag- 
ner and Scidman^ estimate the migration enthalpy to 
be between 0.24 and 0.32 eV, with He being immobile up 
to temperatures of at least 90K, which is consistent with 
the value of 0.28 eV found for 3 He migration by Amano 
and Seidman 70 . The discrepancy between ab initio and 
experiment was explained by Becquart and Domain as 
being due to the presence of strong He-He binding, as 
found in their ab initio calculations, resulting in the for- 
mation of less mobile interstitial He clusters for all but 
the lowest concentrations^ 7 -. This is consistent with the 
work of Soltan et alQ, who found He to be mobile in W 
and Au at temperatures below 5K but with increasing 
suppression of mobility as a function of He concentra- 
tion. 

To this data we add the results of our own investigation 
into He migration in Ni. Following on from the results in 
afmD and afml Fe, we make the reasonable assumption 
that the most stable off-centre octa He configuration is 
a good candidate for the transition state for interstitial 
He migration. The additional uncertainty on the inferred 
migration barrier height from this assumption should be 
less than 0.01 eV. From the results presented in Table ITO1 
this is the off-centre octa {110} configuration, with a cor- 
responding migration barrier height of 0.13 eV, which 
compares well with the experimental value of 0.14 ± 0.03 
eV measured by Philipps and SonnenbergZ^, correspond- 
ing to a migration activation temperature of 55 ± 10 K. 
This barrier height also compares well with our best es- 
timate for austenitic Fe. We therefore tentatively sug- 
gest that the barrier height for interstitial He migration 
in austenitic Fe-Ni based alloys will be in the range 0.1 
to 0.2 eV, resulting in free, three-dimensional diffusion 
well below room temperature. We accept that there is a 
very real possibility of significant local composition de- 
pendence in these concentrated alloys but we speculate 
that the effective barrier height will still be in the given 
range. 



2. Substitutional He migration 

The diffusion of substitutional He generally proceeds 
via the dissociative and vacancy mechanisms 7 - ' 1 3 1 73 , Di- 
rect exchange with a neighbouring solute atom provides 
an alternative mechanism 7 ^, but is highly unlikely to con- 
tribute significantly to diffusion due to the large acti- 
vation energy for the process. For example, our best 
estimate of the barrier height in Ni is 3.50 eV, which 
compares well with that found using an embedded atom 
model (EAM) potential of 3.1 eV by Adams and Wolfer 73 
and means that substitutional He is, essentially, immo- 
bile. 

For many applications, substitutional He is best con- 
sidered as an interstitial He atom strongly bound to a 
vacancy point defect, with a binding energy, ^(He 1 , V). 



The dissociative mechanism for substitutional He mi- 
gration proceeds by the dissociation of He from a va- 
cancy followed by interstitial migration until it becomes 
trapped in another vacancy. As such, the diffusion coef- 
ficient by this mechanism is inversely proportional to the 
vacancy concentratio n 7 ' 8 ' 73 . If thermal vacancies domi- 
nate, the activation energy is given by 7 - ^ 13 ' 73 

Ef ss - = ^(He 1 ) + £ b (He\ V) - E { (V), (5) 

where £' m (He I ) is the migration energy for interstitial He. 
However, if there is a supersaturation of vacancies, for ex- 
ample under irradiation, then the diffusion is dominated 
by the dissociation step and 

Ei iss - = E m (He 1 )+E h (Re 1 ,V), (6) 

which is, essentially, the dissociation energy for substi- 
tutional He from its vacancy, and the diffusion coeffi- 
cient will remain inversely proportional to the vacancy 
concentration^ 3 -. 

The diffusion of a substitutional solute by the vacancy 
mechanism in an fee lattice is usually well described by 
the 5-frequency model of Lidiard and LeClair o 74 ' 75 . A 
key assumption of this model is that when a vacancy 
binds at Inn to a substitutional solute, the solute remains 
on-lattice. However, this is not the case for He, which 
we find relaxes to a position mid-way between the two 
lattice sites to form a V 2 He complex. The possibility 
of solute-vacancy exchange at 2nn is also not included 
in this model, a point we will return to in the following 
discussion. 

Given the strong binding between a vacancy and sub- 
stitutional He at Inn, which we discuss in Section IIV Al 
we assume that the migration of the V2He complex, as 
a single entity, dominates the diffusion by the vacancy 
mechanism a , with a migration energy, i? m (V2He). The 
diffusion coefficient will be proportional to the V2He con- 
centration, which is in turn proportional to the vacancy 
concentration and depends on the binding energy be- 
tween a substitutional He and a vacancy, i?b(He s 7 V). 
The resultant activation energy for substitutional He mi- 
gration by the vacancy mechanism is given by 7 -, 

ET C - = £ m (V 2 He) - £ b (He s , V) + E f (V) (7) 

when thermal vacancies dominate and by 

ET C - = E m (V 2 Hc) - E h (He s , V) (8) 

when there is a supersaturation of vacancies^. 

We have determined the migration energies for the 
V2He complex using a combination of NEB and single 
configuration calculations. In afmD and afml Fe, where 



a Other processes may contribute and a model, generalising the 
five-frequency model of Lidiard and LeClaire, is needed for this 
situation. 
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FIG. 4: Configurations for A-B pairs of interacting substitu- 
tionally sited solutes and defects in afmD Fe. Species A is 
shown in black and species B in grey along with the config- 
uration label. Fe atoms are shown in white with arrows to 
indicate the local moments. Coordinate axes are as in Fig. [1] 

more than one distinct V2He complex exists, we have 
calculated the migration energy along all of the distinct 
paths where the migrating Fe atom retains the sign of 
its magnetic moment. In previous work 4 ^, we found un- 
realistically high migration barriers along paths where 
the moment changed. We label the migration paths for 
V2He migration by the initial and final configurations, 
which are defined in Fig. [4j and present the correspond- 
ing migration energies in Table IVIIII 



Path afmD Fe afml Fe Ni_ 

lb-s-lb 1.033 as lc-S-lc 1.197 

lc-Hc 0.910 0.898 

la->lb 1.216 
lb->la 1.211 



TABLE VIII: Migration energies, _B m (V 2 He), in eV, for the 
V2He complex. The migration paths are labelled by the initial 
and final configurations, as defined in Fig. [4] 

The migration energies lie approximately 0.2 eV higher 
than those for the corresponding single vacancy migra- 
tion in afmD and afml Fe 4 ^ and in Ni, where we found 
a vacancy migration energy of 1.06 eV, in good agree- 
ment with other DFT calculation s 22 ' 47 and with the ex- 
perimental average valued of 1.04 ± 0.04. We suggest 
that this results from the additional energy required to 



move the He atom from its central position in the V2He 
complex back towards the lattice site during migration, 
as observed in all cases. We also contrast these results 
with those for divacancy migration. In afmD Fe, afml 
Fe and Ni we find migration energies for the divacancy 
along the lb->lb path of 0.370, 0.221 and 0.473 eV, re- 
spectively, which are significantly lower than those for 
the V2He complex. In this case the difference arises not 
only from the energy required to move He to an on-lattice 
site during migration but also from its hindrance of the 
migrating Fe atom. 

Vacancy-He exchange at 2nn provides an alternative 
migration path for substitutional He to that of V2He 
migration. We found energy barriers for 2nn exchange 
as low as 0.47 and 0.55 eV in afmD and afml Fe, re- 
spectively, and a value of 0.94 eV in Ni. While these 
results are lower than the migration energies for V2He, 
the repulsive interactions between a vacancy and sub- 
stitutional He at 2nn [see Section IIV A| mean that the 
equilibrium concentrations of such configurations will be 
significantly lower than the V2He concentration, result- 
ing, we believe, in a much lower contribution to total 
diffusion. While this does strengthen our position that 
V2He migration dominates substitutional He diffusion by 
the vacancy mechanism, a model including all the rele- 
vant migration paths is necessary to answer this question 
conclusively. 

Using the results presented here and in Section IIV1 we 
evaluate the expressions in Eqns. <j5j) , ([6]) , (JTJ) and ((5J and 
present the results in Table lIXl 





afmD Fe 


afml Fe 


Ni 


Ef ss -, Eq. © 


0.599-0.788 


0.853-0.902 


1.405 


£l iss ', Eq. © 


2.411-2.600 


2.810-2.859 


2.756 


£X ac \ Eq. © 


2.066-2.413 


2.232-2.251 


2.192 


Eg"; Eq. © 


0.254-0.601 


0.275-0.294 


0.841 



TABLE IX: Activation energies for substitutional He migra- 
tion, in eV, by the dissociative, E^ BB ' , and vacancy, Eg c ; 
mechanisms for thermal (Eqns. ((5j and Q) and supersatu- 
rated (Eqns. © and ©) vacancy concentrations. For afmD 
and afml Fe we give the range of possible values corresponding 
to the distinct migration paths in these states. 

The results clearly differentiate between the two mech- 
anisms and show a strong correlation to corresponding 
results in bcc Fe 25 . When thermal vacancies dominate 
we predict that diffusion will proceed predominantly by 
the dissociative mechanism. If a supersaturation of va- 
cancies exists then the vacancy mechanism clearly has 
the lowest activation energy. However, the vacancy con- 
centration also plays a critical role in determining which 
mechanism dominates through the distinct way it enters 
the expressions for the diffusion coefficients. For suffi- 
ciently low but still supersaturated vacancy concentra- 
tions the dissociative mechanism may become dominant. 
This is, however, most likely to be the case at low tern- 
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peratures where diffusion by either mechanism is likely to 
be negligible. As such, we suggest that vacancy mediated 
diffusion is the most important mechanism in conditions 
of vacancy supersaturation. 

For the case of Ni, Philipps and Sonnenberg& find 
an activation energy for He diffusion of 0.81 ± 0.04 
eV from isothermal, He-desorption spectrometry experi- 
ments. They attribute this result to the diffusion of sub- 
stitutional He by the dissociative mechanism, hindered 
by thermal vacancies, from which they infer an energy 
for dissociation of He of 2.4 eV. Our results agree that 
substitutional He migration will proceed by the dissocia- 
tive mechanism in a thermal vacancy population. There 
is, however, a 0.6 eV difference between our calculated 
activation energy [Table IIX| and experiment. We also 
find a dissociation energy for He from the substitutional 
site of 2.756 eV, which is in excess of the inferred exper- 
imental value. This large discrepancy suggests that the 
inferred experimental mechanism may not be correct. Ab 
initio calculations show that interstitial He atoms bind 
strongly to one another in Ni22. As discussed earlier, just 
such a mechanism was responsible for the suppression of 
interstitial He migration in W and may also explain the 
experimental result in Ni. 



3. Interstitial C and N migration 

The migration of interstitial C and N in both Fe and 
Ni goes from octa site to adjacent octa site. In afmD 
Fe, there are three distinct migration paths, depending 
on where the initial and final octa sites lie. We label 
these as "in-plane", when the octa sites lie in the same 
magnetic plane, "uu" , when the octa sites lie in adjacent 
magnetic planes with the same sign of magnetic moment 
and "ud" , when the octa sites lie in adjacent magnetic 
planes with the opposite sign of magnetic moment. In 
afml Fe, only the "in plane" and "ud" paths are dis- 
tinct and in Ni, all paths are equivalent. Each of these, 
distinct, migration paths will be symmetrical about an 
intermediate state lying in the plane that bisects the di- 
rect path between the two octa sites. In what follows, 
we consider the tetra and (110) crowdion sites as can- 
didate intermediate states. Possible migration paths for 
in-plane migration in afmD Fe are shown in Fig. [5j as an 
example. 

For C, the results in Tables U and [H] show that the 
crowdion configurations are the lowest lying of the pos- 
sible intermediate states. We have performed NEB cal- 
culations in afmD Fe for C migrating from the octa site 
to all of the distinct crowdion sites in order to determine 
the energy profiles along these paths. We find a single 
maximum in the energy at the crowdion configurations. 
We find this is also the case for N migration via the crow- 
dion configuration in Ni, as will be discussed below. On 
this basis and given the significant local dilatation nec- 
essary to form a crowdion, we make the assumption that 
there will be a single energy maximum at all (110) crow- 




FIG. 5: Possible migration paths for interstitial C and N in 
afmD Fe. Paths are shows for migration from initial to final 
octa sites (black circles) lying in the same magnetic plane 
via tetra ud, [110] crowdion and tetra uu intermediate sites 
(grey circles). Fe atoms (white circles) are shown with arrows 
to indicate the local moments. Migration between octa sites 
in adjacent magnetic planes have not been show for clarity. 
Coordinate axes are as in Fig. [T] 



dion sites so that the MEPs and barrier heights for C 
migration in afmD and afml Fe and in Ni can be deter- 
mined from the data in Tables Q] and UU The same can 
also be said for N migration in afmD Fe along uu and 
ud paths. For all other cases, however, the tetra sites 
are lower in energy and we have performed NEB calcula- 
tions with climbing imaged to investigate the migration 
energy profiles along these paths. 

In afmD Fe, our calculations confirm that the tetra 
ud site is the energy barrier for N migration. In afml 
Fe, however, there is evidence of a shallow minimum, 
0.015 eV deep, around the tetra configuration. Results 
in Ni, by contrast, show a clear double-peaked structure 
in the energy profile. We present the results in Fig. [6] 
and include the results for migration via the crowdion 
site for comparison. The results show that the tetra N 
configuration is a stable local minimum, with a depth of 
0.273 eV relative to the transition state, and not a saddle 
point, like the crowdion configuration. Despite this, the 
MEP for N migration is still via the tetra site. 

We summarise our results for the energy barriers and 
MEPs for interstitial C and N migration in Table [X] In 
the Fe reference states there is a significant spread in 
the migration barrier heights for C migration, both along 
distinct migration paths and between the two states. In- 
plane migration clearly exhibits a higher energy barrier 
in both states, which results directly from the tetrago- 
nal distortion of the lattice and the subsequently higher 
energy necessary to form the [110] crowdion transition 
state. The data also suggests a significant dependence on 
the local magnetic order, just as was seen for interstitial 
He migration. The large spread in barrier heights means 
we cannot make any definitive predictions, except that 
diffusion is three-dimensional. However, in any thermo- 
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C 


migration 










N migration 




Path 




afmD Fe 




afml Fe 




Ni 




afmD Fe 


afml Fe 


Ni 






TS/I 


Em 


TS/I 


Em 


TS/I 


Em 


TS/I 


Em TS/I 


Em TS/I 


in plane 


2.033 


[110] crow. 


2.445 


[110] crow. 


1.628 


(110) crow. 


1.558 


tetra ud 


1.899 tetra* 


1.296 tetra* 


uu 


1.443 


[Oil] crow, uu 




as ud 


as 


in plane 


1.602 


[011] crow. 


uu as ud 


as in plane 


ud 


1.310 


[Oil] crow, ud 


2.113 


[Oil] crow, ud 


as 


in plane 


1.384 


[011] crow. 


ud 1.899 tetra* 


as in plane 



TABLE X: Migration energy barrier heights, Em, in eV for interstitial C and N migration in afmD and afml Fe and in Ni. 
Migration is between adjacent octa sites via a transition state/intermediate (TS/I), which is specified in the table, along all of 
the distinct paths for each particular reference state. Where the transition state/intermediate is only an intermediate state on 
the migration path, its name has been marked with an asterisk. 




Reaction Coordinate 

FIG. 6: Formation energy difference, AEj, to the lowest en- 
ergy octa configuration for N migration in Ni via tetra and 
(110) crowdion intermediate states. 

dynamic average, the lower-energy paths will dominate, 
which suggests an effective barrier height around 1.4 eV 
in afmD Fe and 2.1 eV in afml Fe. The afmD Fe value 
is reasonably consistent with the experimentally deter- 
mined activation energies for C migration in austenite of 
1.626 eV-S and 1.531 eV-S. 

In Ni, we find that C migrates via the crowdion site 
with an energy barrier height of 1.63 eV. This con- 
trasts with the 32 atom cell, where the tetra pathway 
is preferred^. Once again, this demonstrates the inad- 
equacy of using a 32 atom cell for solute migration in 
fee Fe and Ni. Experimental results, using a variety of 
techniques applied both above and below the Curie tem- 
perature, Tc = 627 K, for Ni, yield activation energies in 
the range 1.43 to 1.75 eV— , consistent with our results. 
The experimental results also suggest that the influence 
of magnetism on the migration barrier (and enthalpy of 
solution) for C is no more than 0.2 eV and suggests this 
is the likely error in using fm Ni results to estimate those 
in the paramagnetic state. 

Experimental results for C in Fe-Ni austenitic alloys, as 
discussed by Thibaux et alJ^, show only slight changes in 



C mobility as a function of Ni composition in the range 
from 20 to 100 wt% Ni. They also report an activa- 
tion energy of 1.30 eV in an Fe-31wt%Ni austenitic alloy. 
Overall, our results, in conjunction with the experimen- 
tal results we have discussed, suggest that the migra- 
tion energy barrier for C migration will lie in the range 
1.5±0.2 eV across the whole composition range for Fe-Ni 
austenitic alloys. 

For N, the migration barrier lies between 1.38 and 1.60 
eV in afmD Fe, with a value of 1.90 eV in afml Fe. As 
with C, the afmD Fe results are, on average, lower than 
those for the afml state. The result of an Arrhenius fit to 
combined experimental diffusion data for N migration in 
austenite gave a similar value of 1.75 eV— . In Ni, we find 
a barrier height of 1.30 eV, which is in excess of the ex- 
perimental activation energy reported by Lappalainnen 
and Anttila^ of 0.99 ± 0.12 eV. In light of the significant 
variation in experimental results for C migration in Ni, 
these two results are in reasonable agreement and cer- 
tainly to within the 0.2 eV we have suggested earlier as 
a likely error when using ferromagnetic Ni to model the 
paramagnetic state. Overall, these results show that N 
migrates with a significantly lower barrier in Ni than in 
austenitic Fe and we would expect to find an intermediate 
value in Fe-Ni based alloys, more generally. 



C. Solute-solute interactions 

We have performed calculations to investigate the in- 
teractions between pairs of He atoms in substitutional 
and tetra sites in afmD and afml Fe. Configurations 
with single substitutional and tetra-sited He atoms at up 
to 2nn separation were found to consistently relax to a 
vacancy containing two He atoms. While this does not 
yield any useful binding energy data it does indicate that 
there is little or no barrier for this process and places a 
lower limit on the capture radius of a substitutional He 
of around 3 A. Results for pairs of interacting substitu- 
tional and tetra-sited He atoms, as identified in Figs. |31 
and [7j respectively, are given in Table IXII 

Substitutional He pairs show consistent results across 
the Fe reference states with a strong positive binding at 
Inn and slightly repulsive interactions at 2nn [Table [XTj . 
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FIG. 7: Configurations for interacting tetra-sited solutes in 
afmD Fe. Configurations are labelled by the indices attached 
to the appropriate solute atom positions, shown in black. Fe 
atoms are shown in white with arrows to indicate the local 
moments. Coordinate axes are as in Fig. [T] 





afmD Fe 


afml Fe 


A-B/Config. 


Bt 


E h 


Ei E h 


S-S/la 


7.115 


0.934 


7.423 0.946 


S-S/lb 


7.112 


0.937 


as S-S/lc 


S-S/lc 


6.976 


1.073 


7.419 0.950 


S-S/2a 


8.197 


-0.149 


8.570 -0.201 


S-S/2b 


8.109 


-0.060 


8.493 -0.123 


T-T/l-2 


8.614 


0.313 


9.692 0.242 


T-T/l-3 


8.831 


0.096 


as T-T/2-4 


T-T/2-4 


8.215 


0.712 


9.463 0.472 


T-T/l-4 


rlx T-T/2-5 


9.403 0.531 


T-T/l-5 


8.700 


0.227 


as T-T/2-6 


T-T/2-6 


8.643 


0.284 


9.428 0.506 


T-T/2-5 


8.487 


0.440 


9.340 0.594 



TABLE XI: Formation and binding energies in eV for inter- 
acting pairs of He atoms in substitutional (S) and tetra (T) 
sites in afmD and afml Fe. The configurations are labelled 
as in Figs. 0] and [7] for S-S and T-T pairs, respectively. In 
the afmD reference state the binding energies between tetra- 
sited pairs of He atoms have been calculated relative to two 
isolated tetra ud He. For interacting pairs of tetra-sited He 
atoms at 2nn and 3nn separation the configurations are la- 
belled by the initial He positions, which due to the significant 
displacements under relaxation should not be taken as the fi- 
nal positions. Eshelby corrections for S-S pairs were found to 
be negligible but were -0.09 eV for T-T pairs with a resulting 
increase in the T-T binding energies of up to 0.05 eV. 



In our calculations, He atoms at Inn relax directly to- 
wards one another by between 0.38 and 0.44 A, resulting 
in a consistent He-He separation of between 1.67 and 
1.69 A. While still close to the lattice sites, these dis- 
placements are in stark contrast to the insignificant dis- 
placements observed at 2nn. Substitutional He pairs in 



Ni are similar: At Inn the He atoms are displaced to- 
wards one another by 0.37 A to a He-He separation of 
1.75 A. The resultant binding energy, at 0.657 eV, is less 
than in Fe but is in similar proportion to the substitu- 
tional He to vacancy binding energy [see Section IIV A| . 
At 2nn He remains on-lattice with a repulsive binding 
energy of -0.16 eV. 

Pairs of tetra-sited He atoms exhibit significant bind- 
ing energies of up to 0.7 eV in afmD Fe and 0.6 eV in afml 
Fe. Such strong interactions are consistently observed in 
bcc and fee metals. Previous ab initio calculations found 
binding energies of 0.47 eV in Ni22, 0.7 eV in Pd^i and 
A122, 0.4 eV in bcc Fe2£ and 1.0 eV in W 27 . At Inn sepa- 
ration, the He atoms in afmD and afml Fe are displaced 
from the tetra sites only slightly under relaxation. The 
resulting He-He "bonds" all lie along one of the axes of 
the unit cell with lengths in a small range from 1.62 to 
1.68 A, which is consistent with those found for substi- 
tutional He pairs at Inn. At 2nn and 3nn the He atoms 
displace significantly towards one another under relax- 
ation from the initial tetra sites, resulting in He-He sep- 
arations from 1.51 to 1.65 A. These displacements are 
sufficiently large to take the He atoms either to the edge 
of their initial tetrahedral regions or into the adjacent 
octahedral region via one of the faces of the tetrahedron. 
This is most pronounced for the 3nn T-T/2-5 configu- 
ration, which in afml Fe relaxed to a configuration with 
the He atoms within the octahedral region and symmet- 
rically opposite the central position along the [111] axis. 
The situation is similar for the afmD state but one He 
atom is significantly closer to the central position. It 
is worth noting that this is the most stable configura- 
tion in afml Fe and the second most stable in afmD Fe. 
The large binding energies result, simply, from the co- 
operative dilatation of the lattice and the reduction of 
repulsive He-Fe interactions, which are naturally great- 
est when the two He atoms are in close proximity. The 
results at 2nn and 3nn separations show that the local 
dilatation of the lattice around a single interstitial He re- 
sults in an attractive force to other interstitial He atoms 
up to at least 3 A away and encourages the formation of 
clusters. 

To investigate interstitial cluster formation further we 
have determined the most stable configurations with 
three and four He atoms in afmD and afml Fe. For a fixed 
number of He atoms we found many distinct configura- 
tions with similar energies but the most stable clusters 
were reasonably predictable from a simple pair interac- 
tion model, given the data in Table IXII In afmD Fe, the 
most stable He3 configuration found had two He atoms 
in a 2-4 formation [see Fig. [7] with the third occupying 
the nearest octa site. In afml Fe, the most stable was an 
L-shaped 1-2-3 cluster. In the most stable He4 clusters, 
all He atoms occupied tetra sites in a rectangular-planar 
formation with Inn edges, such as a 1-2-3-4 cluster. This 
is, in fact, the most stable arrangement found in afml Fe, 
whereas in afmD Fe a square-planar configuration with 
all He atoms in tetra ud sites was the most stable. The 
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total binding energies for the most stable clusters are 
given in Table IXTlI along with results in Ni 22 . 



Cluster 


afmD Fe 


afml Fe 


Ni 


He 2 


0.712 


0.594 


0.47 


He 3 


1.537 


1.374 


1.25 


He 4 


2.637 


2.561 





TABLE XII: Total binding energies, in eV, for the most stable 
interstitial He clusters containing up to 4 He atoms. Results 
in Ni are from the work of Domain and Becquar1j2S. Eshelby 
corrections were found to be -0.19 and -0.34 eV for He3 and 
He4 clusters, respectively, with corresponding increases in the 
binding energies of 0.14 and 0.27 eV. 

The strong clustering tendency of interstitial He is 
clearly demonstrated by the data. Application of the 
Eshelby corrections only enhances this effect. The bind- 
ing energy for an additional He, that is, £b(He n ) — 
I?b(He n _i), increases with n for the small clusters stud- 
ied here. We would expect this, however, to plateau to an 
additional binding energy of around 1 eV per He atom in 
afmD and afml Fe and in Ni, given that the co-operative 
dilatation of the lattice that gives rise to the binding 
happens locally. Such strong clustering can, not only, re- 
sult in an effective reduction in interstitial He mobility as 
He concentration increases but is also a critical first step 
in the spontaneous formation of Frenkel-pair defects, as 
observed in gold^. Indeed, the most stable He4 configu- 
rations found here show a significant displacement of the 
nearest Fe atom to the cluster off-lattice by 0.94 A in 
afmD Fe and 1.36 A in afml Fe. We consider this pos- 
sibility further in Section TlVBI in the context of V m He„ 
clustering. 



afmD Fe afml Fe 



A-B/Config. 


E { 


E h 


E { 


E b 


C-C/la 


-17.490 


-0.104 


-17.572 


-0.141 


C-C/lb 


-17.487 


-0.106 


N/A 




C-C/lc 


-17.559 


-0.034 


-17.561 


-0.151 


C-C/2a 


-17.420 


-0.174 


-17.487 


-0.226 


C-C/2b 


-17.617 


0.023 


-17.655 


-0.058 


N-N/la 


-17.010 


-0.195 


-17.037 


-0.204 


N-N/lb 


-17.035 


-0.170 


N/A 




N-N/lc 


-17.131 


-0.074 


-17.020 


-0.221 


N-N/2a 


-17.054 


-0.150 


-17.075 


-0.167 


N-N/2b 


-17.212 


0.008 


-17.172 


-0.070 


TABLE XIII: 


Formation and binding 


energies, in eV, 


for in- 



teracting pairs of octa-sited C and N interstitials. The con- 
figurations are as labelled in Fig. [8] Eshelby corrections were 
-0.03 and -0.06 eV in afmD and afml Fe, respectively, with 
resultant increases in the binding energy of 0.02 and 0.03 eV. 

Interactions between pairs of octa-sited C and N atoms 




FIG. 8: Configurations for A-B pairs of interacting octa-sited 
interstitials in afmD Fe. Species A is shown in black and 
species B in grey along with the configuration labels. Fe atoms 
are shown in white with arrows to indicate the local moments. 
The lowest symmetry afmD state is shown to uniquely identify 
all of the distinct configurations. Some of these configurations 
will be symmetry equivalent in the afml state. Coordinate 
axes are as in Fig. \T\ 

at up to 2nn separation in afmD and afml Fe are given in 
Table IXIIII The interactions are, generally, repulsive at 
both Inn and 2nn, with a reasonable consistency between 
the two reference states, although the repulsion is slightly 
stronger in the afml state. For C, the pair binding energy 
is between -0.1 and -0.15 eV at Inn and more repulsive 
at 2nn at up to around -0.2 eV. By contrast, N pairs 
exhibit stronger repulsion than C at Inn, at around - 
0.2 eV and weaker repulsion at 2nn, where the binding 
energy is at most around -0.15 eV. Calculations for C-C 
pairs at up to 4nn separation in afml Fe found a maximal 
binding energy 0.02 eV. We conclude that there will be 
no appreciable positive binding of C-C and N-N pairs in 
our reference states for austenite. 

Experimental determinations of C-C and N-N inter- 
action energies in austenite arc discussed in a review 
by Bhadeshia^. If quasichemical theory, which only in- 
cludes Inn interactions, is used to interpret the existing 
thermodynamic data then a pair binding energy of -0.09 
eV is found for C and -0.04 eV for N. Our results for C in 
afmD and afml Fe are in good agreement with this value 
and while we do find a repulsive interaction between N-N 
pairs, we find a stronger repulsion than for C, which is 
contrary to the results of this analysis. A more detailed 
analysis can be performed using Mossbauer spectroscopy 
data to study the distribution of C atoms in the Fe ma- 
trix, which can be compared with the results of Monte 
Carlo simulations to determine the solute interaction en- 
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ergies at Inn and 2nn 85 . Such an analysis yields C-C 
binding energies of -0.04 eV at Inn and less than -0.08 
eV at 2nn and N-N binding energies of -0.08 eV and -0.01 
eV at Inn and 2nn, respectively^. Our results follow the 
same pattern for the relative strengths of repulsion but 
are significantly in excess of the results of this analysis. 
The agreement is still impressive, however, given the level 
of extrapolation between our two ordered magnetic states 
at OK and temperatures where paramagnetic austenite is 
stable. 

For C in fm Ni, Siegel and Hamilton 68 found C-C bind- 
ing energies at Inn and 2nn of 0.01 and -0.01 eV, re- 
spectively, using comparative DFT calculations to those 
performed in this work. They, furthermore, show that 
this negligible level of binding is consistent with the ex- 
perimental estimates of the C-C pair concentration as a 
function of total C concentration 8 - 6 -. 

From the data presented above we would suggest that 
C-C and N-N interactions in Fe-based austenitic alloys 
will be repulsive at Inn and 2nn, with binding energies 
in the range from -0.1 to -0.2 eV. We would, furthermore, 
expect the level of repulsion to be reduced as a function 
of increasing Ni concentration. 

D. Interactions with substitutional Ni and Cr 
solutes in Fe 

As an initial step in the investigation of the interactions 
of He, C and N with substitutional Ni and Cr solutes in 
austenite we have calculated the formation energies for 
single substitutional Ni and Cr and present the results in 
Table IXTVl b On this basis, the results of our calculations 
of the interactions between He, C and N solutes and sub- 
stitutional Ni and Cr solutes in afmD and afml Fe are 
presented in Table IXVl 



Config. 


fct afmD 


fct afml 




Et /i 


E t fi 


Sub. Ni 


0.083 0.039 


0.145 -0.301 


Sub. Cr 


0.263 0.843 


0.061 1.120 



TABLE XIV: Formation energies, Et, in eV and magnetic 
moments, ji, in \ib for substitutional Ni and Cr solutes in 
austenitic Fe. The sign of the moments indicates whether 
there is alignment (positive) or anti-alignment (negative) with 
the moments of the atoms in the same magnetic plane. 

In both Fe reference states, substitutional He binds 
weakly to Ni, by around 0.1 eV, and has a repulsive in- 
teraction with Cr of -0.1 eV. The similarity to vacancy- 
substitutional Ni/Cr binding is striking^. The similarity 



b These results agree to within errors with our previous work2£ at 
a lower plane-wave cutoff energy and direct the reader there for 
further discussion. 
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TABLE XV: Formation and binding energies, in eV, for sub- 
stitutional Ni/Cr (species A) to substitutional He, tetra He 
and octa C/N (species B) with configurations labelled as in 
Figs . l4l l9l and [TOl respectively. Eshelby corrections to E{ were 
found to be -0.02 eV when interstitial solutes were present but 
negligible for all other quantities. 



also extends to the local moments on Inn atoms sur- 
rounding the substitutional He and vacancy, as was dis- 
cussed in Section Mil These results are also consistent 
with Ni and Cr acting as slightly oversized and undersized 
solutes, respectively, when interacting with point defects 
in afmD and afml Fe, as discussed previously^. We 
would expect the interactions of other transition metal 
solutes with substitutional He to be readily inferred from 
their interactions with vacancies. 

Interstitial He binds weakly to Ni by, on average, 0.09 
eV at Inn and 0.05 eV at 2nn in the Fe reference states. 
We also observe weak positive binding with Cr, but only 
at Inn, where the binding energy is, on average, 0.13 eV. 
Closer observations of the configurations revealed that He 
relaxed slightly away from Ni, but toward Cr at Inn. Ni 
also remained closer to the lattice site than Cr. These ge- 
ometrical results are consistent with Ni and Cr behaving 
as oversized and undersized solutes, respectively, despite 
both exhibiting binding to interstitial He, although the 
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FIG. 9: Configurations for interactions between a substitu- 
tionally sited species (A) and a tetrahedrally sited species 
(B) in afmD Fe. Species A is shown in black, and species 
B in grey. Configurations are labelled by the position of the 
tetrahedrally sited species, as shown. Fe atoms are shown in 
white with arrows to indicate the local moments. Coordinate 
axes are as in Fig. Q] 



binding to Cr is marginally greater. The level of binding 
suggests that Ni and Cr may act as weak traps for mi- 
grating interstitial He at low concentrations. However, 
with increasing concentration and, therefore, likelihood 
that He remains in similar local environments as it mi- 
grates, a direct study of the local composition depen- 
dence of the migration energy becomes necessary. From 
the binding energy data we can speculate, however, that 
such a dependence will also be weak and maintain our 
earlier suggestion that the activation energy for intersti- 
tial He migration will lie in the 0.1 to 0.2 eV range in 
concentrated Fe-Cr-Ni alloys. 

The interactions of octa C and N with substitutional Ni 
and Cr are reasonably consistent in both afmD and afml 
Fe. For C, interactions with Ni are minimal, although 
slightly repulsive, at Inn, whereas positive binding is ob- 
served with Cr on the order of 0.1 eV. The interactions of 
N are similar to those of C but significantly stronger and 
exhibit a repulsion of around -0.1 eV to Ni and attraction 
to Cr of around 0.2 eV. The repulsive interactions with 
Ni are consistent with the lower solubility of C and, par- 
ticularly, N in fee Ni, compared to afmD and afml Fe [see 
Tables HI and HI], and suggests that the interactions are cu- 
mulative. In the case of Cr, such cumulative interactions 
would encourage the formation of Cr-C/N complexes and 
the precipitation of Cr-carbonitrides, as observed exper- 



4b 

o 




4c 



FIG. 10: Configurations for interactions between a substitu- 
tionally sited species A and octa-sited species B in the fec/fet 
afmD reference state. Species A is shown in black and species 
B in grey along with the configuration labels. Fe atoms are 
shown in white with arrows to indicate the local moments. 
The lowest symmetry afmD state is shown to uniquely identify 
all of the distinct configurations. Some of these configurations 
will be symmetry equivalent in the afml state. Coordinate 
axes are as in Fig. [T] 

imentally in non-stabilised austenitic stainless steels^, 
under conditions where these elements are mobile, that 
is at high temperatures or in irradiated environments. 

IV. SOLUTE INTERACTIONS WITH POINT 
DEFECTS 

In this section we consider the interactions of He, C 
and N with a single vacancy (V), in small vacancy-solute 
clusters, V m X„, and with the [001] self- interstitial (SI) 
dumbbell in afmD and afml Fe and in Ni. We present the 
formation (and binding) energies of the underlying and 
most stable defects and defect clusters in Table IXVTl as 
found previously*^. Pairs of vacancies were consistently 
most stable at Inn separation. The most stable tetra- 
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vacancy cluster consists of a tetrahedral arrangement of 
vacancies at Inn to each other. The most stable tri- 
vacancy cluster is formed from this by placing an atom 
near the tetra-vacancy centre. Finally, the most stable 
hexa-vacancy is an octahedral arrangement of vacancies 
with Inn edges. 



Defect 


afmD Fe 


afml Fe 


Ni 


vacancy 


1.812 


1.957 


1.352 


di- vacancy 


3.443 


3.840 


2.688 




(0.181) 


(0.075) 


(0.016) 


tri- vacancy 


4.790 


5.285 






(0.646) 


(0.587) 




tetra-vacancy 


6.479 


7.097 


4.956 




(0.768) 


(0.733) 


(0.451) 


hexa-vacancy 


8.378 


9.210 


6.865 




(2.493) 


(2.534) 


(1.245) 


[001] SI dumbbell 


3.196 


3.647 


4.135 



TABLE XVI: Formation energies, Ef, in eV, for the vacancy, 
the most stable di-, tri-, tetra- and hexa-vacancy clusters, as 
found by Klaver et aZ.— and the [001] SI dumbbell in afmD 
and afml Fe and in Ni. Total binding energies, Eh, in eV, 
are given for the vacancy clusters in brackets below the for- 
mation energies. The results in Fe are consistent with those 
found previously^ 8 . Results in Ni compare well to other DFT 
calculation a 21 ■ 22 i 47 i S8 i 89 _ Eshelby corrections were found to be 
negligible except for the tetra-vacancy in Ni at -0.03 eV, the 
hexa-vacancy at -0.06 -0.03 and -0.05 eV in afmD Fe, afml 
Fe and Ni, respectively, and the dumbbell at -0.05, -0.08 and 
-0.10 eV in afmD Fe, afml Fe and Ni, respectively. The only 
non-negligible effect on binding energies was for the hexa- 
vacancy, where increases of 0.05 eV, 0.03 and 0.03 eV apply 
in afmD Fe, afml Fe and Ni, respectively. 



A. Vacancy - solute interactions 

We present the formation and binding energies for con- 
figurations containing a single vacancy and solute atom, 
at up to 2nn separation, in Table IXVIII 

1. V-He binding 

We observe strong binding of between 0.60 and 0.66 
eV for V-Sub He pairs at Inn in both Fe reference states. 
This is significantly greater than the binding between 
vacancy pairs^ and represents the simplest case of en- 
hanced vacancy binding by He, as we will discuss in what 
follows. We find that He does not remain on-lattice at 
Inn to a vacancy but relaxes to a position, best described 
as, at the centre of a divacancy. With this perspective, 
the V-Sub He binding represents the significant energetic 
preference of He for the greater free volume available at 
the centre of a divacancy over a single vacancy. At 2nn, 
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TABLE XVII: Formation and binding energies in eV for va- 
cancy (species A) to substitutional He, tetra He and octa C/N 
(species B) with configurations labelled as in Figs.l4"Il9"land ll0l 
respectively. Configurations with a single solute atom in the 
substitutional position (Sub.), where stable, are also consid- 
ered as an interstitial solute interacting with a vacancy. In 
afmD Fe, the vacancy-tetra He binding energies were calcu- 
lated relative to tetra ud He. Binding energies between octa 
C and N solutes and a vacancy at 3nn and 4nn separations 
were investigated but did not exceed 0.03 eV. The only non- 
negligible Eshelby corrections found were for the binding en- 
ergies between a vacancy and interstitial solutes at no more 
than 0.02 eV in magnitude. 



the interactions are repulsive, at around -0.1 eV, which 
is slightly greater than that observed between vacancy 
pairs 48 . He remains on-lattice in these configurations, 
which explains the lack of enhanced binding at 2nn sepa- 
ration. The situation in Ni is very similar, where we find 
binding energies of 0.356 eV and -0.127 eV at Inn and 
2nn, respectively. 

Interstitial He binds strongly to a vacancy to form a 
substitutional He configuration. The same is also true 
in Ni, where we find a binding energy of 2.627 eV, in 
good agreement with previous work™ Configurations 
with tetra He at Inn to a vacancy are unstable. At 
2nn, however, we find stable configurations with weak 
repulsive or attractive binding, depending on the con- 
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figuration. The fact that no stable configurations were 
found with tetra He at up to 2nn from a substitutional 
He atom demonstrates that the addition of a single He 
to a vacancy significantly increases the capture radius for 
interstitial He. We expect this effect to increase with the 
subsequent addition of He, given the additional pressure 
and dilatation that would be exerted on the surrounding 
lattice. 



2. V-C and V-N binding 

C binds to a vacancy by up to 0.38 eV at Inn in the Fe 
reference states and weakly at 2nn. This level of binding 
agrees well with previous experimental and theoretical 
work in austenite and austenitic alloys^. We find that V- 
N binding is significantly stronger than for C with bind- 
ing energies in the range from 0.3 to 0.6 eV at Inn and 
around 0.1 eV at 2nn. For both C and N, the substitu- 
tional configuration is strongly disfavoured. As discussed 
in Section IHl) the substitutional C and N configurations 
in afmD Fe were found to be unstable and the configu- 
ration labelled V-C/Sub in Table IXVHl has C in a stable 
position off-lattice by 0.77 A. Overall, these results bear 
a strong similarity to those found in bcc Fe^l, where bind- 
ing energies of 0.47 and 0.71 eV were found for C and N 
at Inn to a vacancy, respectively. 

Results in Ni are broadly similar to those in Fe. We 
find a V-C binding energy of 0.062 eV at Inn and 0.121 
eV at 2nn. V-N binding is, again, stronger, than C, with 
energies of 0.362 eV and 0.165 eV at Inn and 2nn, re- 
spectively. We also find a strong repulsion from the sub- 
stitutional site. We note that the V-C binding at Inn 
seems anomalously low, given the other results but no 
problems were found with this calculation and other test 
calculations found the same stable structure and energy. 

The significant V-C and V-N binding energies suggest 
that the relatively less mobile solutes could act as va- 
cancy traps, much as was found in bcc Fe^ 2 - ' 37 ' 90 . This 
would certainly be the case if dissociation of the complex 
was required before the vacancy could freely migrate but 
the possibility of cooperative migration also exists. In 
the fee lattice there are many possible migration path- 
ways that would avoid the dissociation of this complex, 
including some that would maintain a Inn separation. A 
complete study of these possibilities is beyond the scope 
of this work but preliminary calculations in Ni show that 
the energy barriers for C and N jumps that would main- 
tain a Inn separation to the vacancy are around half the 
value for the isolated solutes at around 0.75 eV. In con- 
trast, vacancy jumps that maintain a Inn separation were 
found to be significantly higher than those for the isolated 
vacancy but jumps from Inn to 2nn separation and back 
exhibited lower or comparable energy barriers. Whilst 
these calculations are preliminary, they do indicate the 
distinct possibility of cooperative vacancy-solute motion 
that would avoid dissociation of the complex. The im- 
plications for an absence of vacancy pinning and for the 



enhanced diffusion of C and N solutes in the presence of 
vacancies in austenitic alloys makes this an interesting 
subject for further study. 



B. Vacancy - solute clustering 

Small Vacancy-He (V m He„) clusters have been found 
to be highly stable both experimentally ^ 1 1 ' 1 2 and using 
DFT technique s 19 ' 22-25 ' 28 in a number of metals and are, 
therefore, critically important as nuclei for void forma- 
tion. Experimental evidence in bcc Fei£ has also shown 
that C can act as a vacancy trap through the formation 
of small, stable V m C„ clusters, which has been confirmed 
in a number of DFT studies 3 ^ - . 3 !. Small V m N„ clusters 
have also been shown to exhibit similar stability 3 ^. 

In this section we present the results of a large num- 
ber of DFT calculations to find the most stable V m X„ 
clusters, where X is He, C or N, in afmD and afml Fe. A 
comprehensive search for the most stable configuration 
was only practicable for the smaller clusters. For larger 
clusters, a number of distinct initial configurations, based 
around the most stable smaller clusters, were investigated 
to improve the likelihood that the most stable arrange- 
ment was found. The total binding energies for the most 
stable configurations can be found in Table IXVIIII 



1. V m He n clusters 

The geometries of the relaxed V m He„ clusters were 
constrained by the tendency to maximise He-He and Hc- 
Fe separations within the available volume and, there- 
fore, minimise the repulsive interactions. In a single va- 
cancy we found that this led to the following structures: 
two He formed dumbbells centred on the vacancy with 
He-He bond lengths around 1.5 A; three He formed a 
near-equilateral triangle with bond lengths of between 
1.6 and 1.7 A; four He formed a near-regular tetrahedron 
with bond lengths between 1.6 and 1.7 A; five He formed 
a near-regular triangular bipyramid with bond lengths 
between 1.6 and 1.8 A and six He formed a near-regular 
octahedron with bond lengths between 1.6 and 1.8 A. In 
clusters with more than one vacancy, a single He atom 
relaxed to a central position. Additional He tended to 
form similar clusters to those seen in a single vacancy 
but now around the centre of the vacancy cluster. The 
tri-vacancy case is interesting because previous DFT cal- 
culations in austenite^ 8 - found that a configuration con- 
sisting of a tetrahedral arrangement of vacancies with one 
Fe atom near the centre of the void, which can be consid- 
ered as the smallest possible stacking fault tetrahedron 
(SFT), was more stable than the planar defect of three 
vacancies with mutual Inn separations. The addition of a 
single He atom was enough to reverse the order of stabil- 
ity with a difference in the total binding energy of 0.8 eV 
in afml Fe, in favour of the planar defect. We suggest that 
this result should readily generalise, with planar defects 
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(a) afmD Fe, He binding 



(b) afml Fe, He binding 




No. Vacancies, m No. Vacancies, m 

(c) afmD Fe, V binding (d) afml Fe, V binding 




12 3 4 

No. Vacancies, m 




2 3 4 

No. Vacancies, m 



(e) afmD Fe, SI binding 



(f) afml Fe, SI binding 



FIG. 11: Binding energies, in eV, for a He atom, V or SI to an existing cluster to form one with the V m He n stoichiometry in 
afmD Fe (figures (a), (c) and (e)) and afml Fe (figures (b), (d) and (f)) Fe. Interstitial He cluster data has been included in 
figures (a) and (b) for completeness. 
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TABLE XVIII: Total binding energies, E h , in eV for the most 
stable V m X n clusters found in afmD and afml Fe, where X 
is He, C or N. Results in Ni are also given for C and N and 
can be found in the literature for H o 19 ' 22 . Eshelby corrections 
to Eh for V m He n clusters were found to be below 0.05 eV in 
magnitude except for VHes and V m He n with m and n equal 
to 3 or 4, which were below 0.1 eV and VHe6, which was 0.2 
eV. For C and N clusters, the corrections were below 0.02 eV 
in magnitude except for those with six vacancies or with four 
or more N atoms, where the corrections were up to 0.1 eV for 
most but were 0.2 eV for VN6 in afml Fe and VN5 in Ni and 
0.3 eV for VN 6 in Ni. 



being more stable than SFTs with sufficient addition of 
He. That said, however, planar defects have been found 48 
to be less stable than three-dimensional proto-voids and 
this situation is unlikely to change with the addition of 
He due to the greater free volume of the latter clusters. 

In Section IIII A II the addition of a single He to a va- 
cancy was found to have very little effect on the local 
magnetism. The addition of He to vacancy clusters was 
generally found to have very little effect on the total mag- 
netic moment of the supercells containing the cluster. 
The only exception was for the single vacancy in afmD 
Fe, although it took the addition of six He atoms to sig- 
nificantly change the magnetic moment. Even in the ab- 
sence of vacancies, a cluster of at least three He atoms 
was necessary to influence the total magnetic moment. 

In Fig. Qj] we present results for the binding energy 
of either a He atom, vacancy (V) or [001] self-interstitial 
dumbbell (SI) to an already existing cluster to form one 
with the V m He n stoichiometry. These results show that 
He consistently binds strongly to an existing cluster and 



that the strength of the binding only increases with m. 
For a fixed value of n, this binding energy will converge to 
the formation energy for interstitial He [see Table H] as m 
increases and is well on the way to doing that for n = 1. 
For fixed m the additional He binding energy appears to 
plateau as n increases although it should diminish even- 
tually as the pressure within the cluster builds. 

The binding energy for an additional vacancy is also, 
consistently, positive. The presence of He significantly in- 
creases this additional binding for all values of m, which 
is consistent with the observation that it aids the nu- 
cleation, stabilisation and growth of voids in irradiated 
environment o 10 i 13 ~ — . For fixed n, the data shows that 
the vacancy binding energy is tending to a plateau as m 
increases and is consistent with the fact that all of these 
curves should converge to the vacancy formation energy. 

The SI binding energy can be related to the vacancy 
binding energy as follows, 

£ b (SI, V m+1 He n ) = £ f (SI) + Ef(V) - E h (V,V m Re n ), 

(9) 

which implies that the spontaneous emission of an SI 
from an existing cluster will be energetically favourable 
if and only if the binding of the newly created vacancy 
is greater than the Frenkel pair formation energy. The 
data shows that the SI binding energy clearly decreases 
as He concentration is increased at fixed m and for suffi- 
ciently high concentration will become negative. Indeed, 
it is energetically favourable for an interstitial He cluster 
with four He atoms in afml Fe, and most likely for five He 
atoms in afmD Fe, to spontaneously emit an SI defect. 
This mechanism was proposed to explain the observa- 
tion of He bubbles in Au samples after sub-threshold He 
implantation 83 and could also explain observations of He 
trapping in Ni2i, where the He was introduced by natu- 
ral tritium decay to avoid implantation-produced defects. 
Our results show that this would, most likely, occur in 
austenite and austenitic alloys and could lead to bubble 
formation, with the potential for blistering in the pres- 
ence of, even low-energy, bombardment by He ions, as 
seen in W 92 ^ 93 . 

As a whole, the binding energy data is qualitatively 
similar to DFT results in AlS 3 . and Pd2^ and is quanti- 
tatively similar to results in bec Fe2£ and N i 19 ' 22 . This 
observation gives us confidence that our results are, not 
only, applicable to austenite but to austenitic alloys more 
generally. 

The binding energy data above has also been used to 
determine the dissociation energy, that is the energy of 
emission, of He, V or SI from a V m He„ cluster using the 
simple ansatz that the dissociation energy, £"di S s.(X), for 
species, X, is given by 

£ diss .(X) = E h (X) + £ m (X), (10) 

where E m (X) is the migration energy for isolated species, 
X. We present results for the dissociation energies in 
Fig. 021 using the migration energies in Table IXTXl 

There is a strong and distinct dependence on the He 
to vacancy ratio, n/m, for the dissociation energies of 
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FIG. 12: Dissociation energies, i?diss.(X), in eV, for species X from a V m He„ cluster, where X is a He, V or SI. Results are 
presented for (a) afmD Fe and (b) afml Fe versus the He to vacancy ratio, n/m. The solid curves are simple polynomial fits to 
the data and are present to aid visualisation. 




Species,X 


afmD Fe 


afml Fe 


He 


0.160 


0.070 


V 


0.743 


0.622 


SI 


0.196 


0.254 



TABLE XIX: Migration energies, _E m (X), in eV, for Species, 
X, where X is He, V or SI. For He, the lowest values from 
Table I VIII were used. For V, the lowest vacancy migration 
energies from Klaver et a/.— were used. The SI migration 
energies were calculated here as that for a [100] dumbbell SI 
migrating between two lattice sites at Inn separation within 
a magnetic plane using identical settings to Klaver et a?.— 



the three species. Both graphs exhibit a clear crossover 
between the He and V curves at around n/m =1.3 and 
another between the He and SI curves at about n/m — 
6. An identical Hc-V crossover ratio was found in bcc 
Fe2£ and fee A122. For n/m below 1.3 the clusters are 
most prone to emission of a vacancy, between 1.3 and 
6 He has the lowest dissociation energy and above 6 SI 
emission is the preferred dissociation product. The slope 
of the curves ensures that emission of the species with the 
lowest dissociation energy will make the resulting cluster 
more stable. At sufficiently high temperatures that these 
processes are not limited by kinetics this should lead to 
the formation of the most stable clusters, which have 
an n/m value at the He-V crossover, where our results 
predict a minimum dissociation energy of around 2.8 eV 
in both afmD and afml Fe. 



2. V m C„ and V m N n clusters 

In fct afmD and afml Fe and in Ni we considered VX„ 
clusters with octa-sited C and N at Inn to the vacancy 
and configurations where C and N are close enough to 
form C-C and N-N bonds within the vacancy. Our re- 
sults for VC2 and VN2 clusters are given in Table IXXI 
We found that VX2 clusters with octa-sited C and N are 
most stable when the C /N atoms are as far apart as pos- 
sible, that is, opposite one another across the vacancy. 
For these configurations the total binding energy is more 
than the sum of the binding energies for each single solute 
to the vacancy, indicating either some chemical or coop- 
erative strain interaction. We found that C-C dumbbells 
centred on the vacancy are stable, with bond lengths be- 
tween 1.38 and 1.48 A, that is, much shorter than the 
separations between octahedral sites. The most stable 
lie along (100) directions and binding over and above 
that for octa sited C was found in afmD Fe and Ni. The 
enhancement in binding upon forming a C-C dumbbell is 
not, however, as pronounced as was seen in bcc Fo 25 ' 34 ~ — . 
We also found stable configurations with N-N dumbbells 
in Fe and Ni with bond lengths between 1.34 and 1.49 A, 
although they exhibit a much lower, and generally neg- 
ative, total binding energy compared to configurations 
with octa sited N atoms. 

For VX3 clusters, we investigated all possible configu- 
rations with three octa-sited C or N solutes in addition to 
those with a C-C dumbbell and an octa-sited C solute in 
one of the four octa sites perpendicular to the dumbbell 
axis and a configuration with three C atoms close enough 
for C-C bonding. Although we do find stable configura- 
tions with C-C bonding in either a dumbbell or triangu- 
lar arrangement in both Fe and Ni, these arrangements 
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afmD Fe 


afml Fc 


Ni 


Config. 


Ei 


E b 


E { 


Eb 


Ef Eb 






VC 2 


clusters 






la-la (opp.) 


-16.226 


0.444 


-16.551 


0.795 


-15.692 0.199 


lb-lc 


-16.166 


0.385 


-16.466 


0.711 


as la-la (opp.) 


la-la (adj.) 


1 ^ Q08 
-Lu. yuo 


n 1 97 


-16.163 


0.408 


-15.396 -0.097 


la-lb 


-15.832 


0.050 
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0.361 


as la-la (adj.) 


la-lc 


-16.104 


0.323 


as la- lb 


as la-la (adj.) 


[100] dumb. 


-16.087 


0.305 


-16.224 


0.469 


-15.915 0.422 


[001] dumb. 


-16.265 


0.484 


-16.262 


0.507 


as [100] dumb. 


[110] dumb. 


-15.407 


-0.375 


-15.238 


-0.517 


-15.380 -0.113 


[Ill] dumb. 


rlx [001 


] dumb. 




-15.565 0.072 






VN 2 


clusters 






la-la (opp.) 


-16.373 


0.981 


-16.579 


1.295 


-14.559 0.872 


lb-lc 


-16.124 


0.732 


-15.932 


0.648 


as la-la (opp.) 


la-la (adj.) 


-16.113 


0.720 


-16.237 


0.953 


-14.279 0.591 


la-lb 


-16.051 


0.658 


-16.037 


0.753 


as la-la (adj.) 


la-lc 


-16.228 


0.835 


as la- lb 


as la-la (adj.) 


[100] dumb. 


-14.451 


-0.942 


-14.487 


-0.796 


-13.616 -0.071 



TABLE XX: Formation and total binding energies, in eV, 
for the interactions of a vacancy with two octa sited C or N 
solutes. Configurations with C or N in octa sites at Inn to the 
vacancy are labelled by the positions of the two solutes as in 
Fig. 1101 When both octa solutes are in the same plane as the 
vacancy the configurations are additionally labelled by their 
relative orientation i.e. opposite (opp.) or adjacent (adj.) to 
one another. Doubly mixed dumbbells centred on the vacancy 
site were also considered as configurations of an interacting 
vacancy with two octa solutes and the total binding energies 
were calculated accordingly. Eshelby corrections to both Ef 
and Eb were found to be no more than 0.03 eV in magnitude. 



are the least stable and exhibit significant, negative total 
binding energies. The most stable arrangements consist 
of three octa sited C atoms placed as far apart as possi- 
ble, for example in three la sites relative to the vacancy 
as in Fig. [TO] However, the total binding energies for the 
most stable VC3 clusters [see Table IXVIIIj are less than 
for VC2, which implies that a vacancy can only bind up 
to two C atoms within a vacancy. A vacancy may still, 
however, bind more than at 2nn octa sites but we did not 
investigate this possibility due to the strongest binding 
being at Inn to the vacancy and due to the large number 
of possible configurations. 

The most stable VN3 clusters have the same geometry 
as found with C but, in contrast, are more stable than 
VN2 clusters. Beyond this point, we found that the total 
binding energy only increases for up to four N atoms in 
afml Fe and Ni but increases all the way up to six N atoms 
in afmD Fe. That said, however, the binding energy per 
N atom only increases up to a VN2 cluster in all reference 
states. The equilibrium concentrations of clusters with 
more than two N atoms, which can be calculated using 
the law of mass-action^, would very likely be negligible, 



even at room temperature. Despite their magnitude, the 
Eshelby corrections do not change these conclusions but 
would result in the total binding energy increasing all the 
way up to six N atoms in afml Fe, as was found for afmD 
Fe. 

We investigated site preference and binding for single 
C and N solutes to the most stable di-, tetra- and hexa- 
vacancy clusters in afmD and afml Fe and in Ni. For the 
V2C cluster in afmD Fe, we considered all Inn octa sites 
to the three distinct types of Inn di- vacancy as well as 
configurations with C at the centre of all Inn and 2nn 
di-vacancy clusters. For the octa sites, C was found to 
bind to existing di-vacancy clusters with similar binding 
energies to a single vacancy, that is with Eb(G, V2) in 
the range from 0.03 to 0.32 eV. The most stable of these, 
which was also found to be the most stable V2C cluster, 
contained the most stable di-vacancy and bound C more 
stably than to a single vacancy. We found that C was 
repelled from the centre of a Inn di-vacancy lying within 
a magnetic plane but bound to the other two Inn di- 
vacancies with similar energies to those found in octa 
sites. As in bcc Fo 3 b 35 ' 36 , the most preferred site for C 
was at the centre of a 2nn di-vacancy, with ^(C, V2) = 
0.35 eV. However, this was not sufficient to overcome the 
difference in stability between Inn and 2nn di-vacancies 
in afmD Fe^ and did not, therefore, form them most 
stable V2C cluster, in contrast to in bcc Fe. 

The analysis above motivated the use of only the most 
stable Inn di-vacancy in the remaining calculations along 
with configurations containing solutes at the centre of 
2nn di-vacancies. For N in afmD Fe, the order of site 
preference mirrors that for C. An N solute is capable of 
stabilising a 2nn di-vacancy configuration but the most 
stable V2N cluster shared the same geometry as for C 
with a binding energy to the underlying di-vacancy of 
0.56 eV, which is, again, in excess of the binding to a 
single vacancy. 

The situation in afml Fe and Ni was found to be rather 
similar to that of afmD Fe. For both C and N, the site at 
the centre of a Inn (in-plane) di-vacancy was disfavoured. 
The most stable configuration generally contained an 
octa-sited solute bound to a Inn di-vacancy. The only 
exception was in afml Fe, where a configuration with C 
at the centre of a 2nn di-vacancy lying within a magnetic 
plane had a greater total binding energy but only by 0.03 
eV. This, most likely resulted from the much smaller en- 
ergy difference between Inn and 2nn di-vacancies in afml 
Fe of compared to afmD Fe^ and to Ni, where we find 
an energy difference of 0.1 eV in favour of the Inn di- 
vacancy. In the most stable clusters, the binding of the 
solutes to the underlying di-vacancy was, once again, in 
excess of the binding to a single vacancy. 

For the binding of C and N to the most stable tetra- 
vacancy, we found that the central position was extremely 
disfavoured. We investigated all configurations with so- 
lute atoms in an octa site at Inn to at least a single 
vacancy. We also performed calculations with solute 
atoms placed initially at random within the proto-void 



but found that these relaxed to octa sites already consid- 
ered. Configurations with only a single vacancy at Inn 
to the solute were found to be the most stable. The total 
binding energies for these configurations are given in Ta- 
ble [STUD Using these results we found that the binding 
of C and N to the tetra-vacancy was in excess of that 
to a di- vacancy and a single vacancy, in all cases except 
for N in afml Fe, where the binding to the tetra-vacancy 
and di-vacancy reversed order, although they only differ 
by 0.02 eV. 

For the hexa-vacancy, the central octa site was un- 
stable for both C and N in afmD Fe. In afml Fe and 
Ni it was stabilised by symmetry but still strongly dis- 
favoured. This repulsion is, however, significantly less 
than was observed for the tetra-vacancy. Closer observa- 
tion showed that while the nearest neighbouring Fe and 
Ni atoms to the solutes moved very little under relaxation 
in the tetra-vacancy, the contraction in bond length was 
between 25 and 30 % in the hexa-vacancy from an ini- 
tial separation of around 3 A. This demonstrates how 
important the formation of strong chemical bonds with 
characteristic bond lengths is to the stability of configu- 
rations containing C and N in Fe and Ni. 

We investigated the stability of configurations with C 
and N in all octa sites at Inn to at least one vacancy in 
the hexa-vacancy cluster. We found that there were ad- 
ditional stable sites, lying along (100) axes projected out 
from the centre of the hexa-vacancy. For C, these sites 
were found to lie between the first vacancy reached along 
these axes and the next octa site out. They are close to 
but distinct from the octa sites and we, therefore, refer 
to them as octa-b sites. For N, stable sites were found 
between the centre of the hexa-vacancy and the first va- 
cancy reached along the (100) axes and we refer to these 
as off-centre sites. We found that C was, consistently, 
most stable in an octa-b site whereas N preferred octa 
sites with two vacancies at Inn, although an off-centre 
site along [001] was the most stable in afmD Fe. 

Once again, the binding energy between the so- 
lutes and hexa-vacancy was greater than for all smaller 
vacancy clusters. We summarise these results for 
E\,(K, V m ) in Fig. [T3l which clearly shows the increase 
in binding energy as the vacancy cluster becomes larger. 
It also clearly shows that in the same reference state, the 
binding energy for N is consistently greater than for C 
and that the binding energies in afml Fe lie above those in 
afmD Fe. The one anomalous point is the binding energy 
for C to a hexa-vacancy in afmD Fe, which is much larger 
than the trends would suggest. Other configurations with 
C in an octa-b site in afmD Fe exhibited similar levels of 
binding and no problems with any of these calculations 
or instabilities in the relaxed structures could be found. 



C. [001] dumbbell SI - solute interactions 

We investigated the binding of He, C and N solutes 
to a [001] dumbbell in afmD and afml Fe and in Ni and 
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FIG. 13: Binding energy, _B b (X,V m ), in eV, where X is C 
(filled symbols and solid lines) or N (unfilled symbols and 
dashed lines) in afmD Fe (red circles), afml Fe (green squares) 
and Ni (blue diamonds). The binding energies were calculated 
for the most stable clusters. 



present the results in Table IXXll 

We found that interstitial He, placed initially Inn to a 
[001] SI dumbbell, either spontaneously displaced under 
relaxation to a 2nn site or exhibited a repulsive binding 
energy in Fe. At 2nn, however, a positive binding en- 
ergy was observed, up to almost 0.2 eV, as was found in 
N i 19 i 22 . Eshelby corrections do not qualitatively change 
these results and would only act to enhance the binding 
at 2nn. This positive binding energy is comparable to 
that in bec Fe2£ but while significant, it is only likely to 
result in mutual trapping at low temperature, given the 
high mobility of the two species. Taken as a model for 
the binding of interstitial He to other overcoordinated 
defect sites, such as near dislocations and grain bound- 
aries, however, this result does show that He would be 
likely to be trapped at such sites, leading to interstitial 
He cluster formation and spontaneous bubble nucleation 
and growth, as discussed earlier. It is worth mentioning 
that bubble nucleation by this mechanism would hap- 
pen much more readily at grain boundaries where, due 
to their disorder, vacancies can be formed without the 
additional SI. 

A substitutional He atom in the lb and lc sites [see 
Fig. 2] to a [001] SI dumbbell resulted in the spontaneous 
recombination of the vacancy and SI and the kick-out of 
an interstitial He atom. At all other Inn and 2nn sites 
except 2b in Ni, however, stable complexes with binding 
energies of up to around 0.2 eV were formed. Barri- 
ers to recombination for these complexes, while positive, 
were not calculated in this work. These results do, how- 
ever, show that substitutional He, and most likely, other 
V m He n clusters can act as trapping sites for SI dumb- 
bells in austenite and austcnitic alloys with a capture 
radius extending out to at least 2nn. We can also spec- 
ulate that, once trapped, recombination will be likely to 
occur. 
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TABLE XXI: Formation and binding energies in eV for [001] 
self-interstitial dumbbell (species A) - solute (species B) in- 
teractions. For octa-sited C and N solutes the configurations 
are labelled as in Fig. 1101 He interactions were investigated 
with He sited substitutional^ and tetrahedrally with config- 
urations labelled as in Figs. [4] and [9] respectively. Configu- 
rations with substitutional He in lb and lc positions relative 
to a [001] dumbbell SI were unstable to defect recombination 
and interstitial He kick-out. In Ni, the [001]-tetra He bind- 
ing energy was observed to be 0.20 eV— . Eshelby corrections 
were found down to -0.1 eV for configurations containing sub- 
stitutional He but the related increases in Eh were no more 
than 0.02 eV. For configurations containing interstitial solutes 
E corr ' could be as low as -0.2 eV with corresponding increases 
in E h up to 0.08 eV. 



Both C and N are either repelled from Inn and 2nn 
sites to a [001] SI dumbbell or show very little positive 
binding, much as was observed in bec Fe^ 4 .. Eshelby cor- 
rections do not change this conclusion in Fe but would 
result in binding of around 0.1 eV at 2nn in Ni. Moti- 
vated by the result that C does exhibit positive binding 
to the most stable SI and small SI clusters in bec Fei2 
at further separation, we investigated this possibility here 
and found sites with binding energies from 0.05 to 0.1 eV 



at 4nn to the dumbbell, which would only be enhanced 
by Eshelby corrections. These sites can be related to the 
corresponding ones in bec Fe by a Bain transformation 94 
and the binding almost certainly results from strain field 
effects in both cases. The fact that such binding was 
found to increase with interstitial cluster sizei^ means 
that Cottrell atmospheres^ of C and N are very likely 
to form around other overcoordinated defects, such as 
dislocations and grain boundaries, in both ferritic and 
austenitic alloys under conditions where these species are 
mobile. 



V. CONCLUSIONS 

An extensive set of first-principles DFT calculations 
have been performed to investigate the behaviour and 
interactions of He, C an N solutes in austenite, dilute Fe- 
Cr-Ni alloys and Ni as model systems for austenitic steel 
alloys. In particular, we have investigated the site stabil- 
ity and migration of single He, C and N solutes, their self- 
interactions, interactions with substitutional Ni and Cr 
solutes and their interactions with point defects typical 
of irradiated environments, paying particular attention 
to the formation of small V m X n clusters. 

Direct comparison with experiment verifies that the 
two-state approach used to model austenite in this work 
is reasonably predictive. Overall, our results demonstrate 
that austenite behaves much like other fee metals and is 
qualitatively similar to Ni in many respects. We also 
observe a strong similarity between the results presented 
here for austenite and those found previously for bec Fe. 

We find that interstitial He is most stable in the tetra- 
hedral site and migrates via off-centre octahedral transi- 
tion states with a migration energy from 0.1 to 0.2 eV 
in austenite and 0.13 eV in Ni. The similarity of these 
results and the weak interactions with Ni and Cr so- 
lutes in austenite suggests a migration energy in Fe-Cr-Ni 
austenitic alloys in the 0.1 to 0.2 eV range. Interstitial 
He will, therefore, migrate rapidly from well below room 
temperature until traps are encountered. Its strong clus- 
tering tendency, with an additional binding energy ap- 
proaching 1 eV per He atom in austenite and 0.7 eV in Ni, 
will lead to a reduction in mobility as interstitial He con- 
centration increases. Interactions with overcoordinated 
defects, which are on the order of a few tenths of 1 eV, 
will result in the build up and clustering of interstitial 
He at dislocations and grain boundaries. The most sta- 
ble traps, however, are vacancy clusters and voids, with 
binding energies of a few eV. The strength of this bind- 
ing means that growing interstitial He clusters eventually 
become unstable to spontaneous Frenkel pair formation, 
resulting in the emission of a self-interstitial and nucle- 
ation of a VHe„ cluster. The binding of additional He 
and vacancies to existing V m He„ clusters increases sig- 
nificantly with cluster size, leading to unbounded growth 
and He bubble formation in the presence of He and va- 
cancy fluxes. The most stable clusters have a helium to 
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vacancy ratio, n/m, of around 1.3, with a dissociation 
energy for the emission of He and V of 2.8 eV in austen- 
ite and Ni. Generally, we assume that V m He n clusters 
are immobile. For the simplest case of substitutional He, 
however, migration is still possible. In a thermal vacancy 
population, diffusion by the dissociative mechanism dom- 
inates, with an activation energy of between 0.6 and 0.9 
eV in Fe and 1.4 in Ni. In irradiated environments, how- 
ever, the vacancy mechanism dominates and diffusion can 
proceed via the formation and migration of the stable 
V2He complex, with an activation energy of between 0.3 
and 0.6 eV in Fe and 0.8 eV in Ni. 

We find that C and N solutes behave similarly, both 
in austenite and Ni, although the interactions of N are 
stronger. The octahedral lattice site is preferred by both 
solutes, leading to a net expansion of the lattice and a 
reduction of the c/a ratio in the afmD and afml Fe ref- 
erence states. Both solutes also stabilise austenite over 
ferrite and favour ferromagnetic over antiferromagnetic 
states in austenite. Carbon migrates via a (110) transi- 
tion state with a migration energy of at least 1.3 eV in 
austenite and of 1.6 eV in Ni. For N, migration proceeds 
via the crowdion or tetrahedral sites, depending on path, 
with a migration energy of at least 1.4 eV in austenite 
and 1.3 eV in Ni. Pairs of solute atoms are repelled at 
Inn and 2nn in austenite and do not interact in Ni. Both 
C and N interact very little with Ni solutes in austenite 
but bind to Cr, which may act as a weak trap and en- 
courage the formation of Cr-carbonitrides in conditions 
where the solutes are mobile. Carbon binds to a vacancy 
by up to 0.4 eV in austenite and 0.1 eV in Ni, with N 
binding more strongly at up to 0.6 eV in austenite and 
0.4 eV in Ni. While this may suggest that C and N act 
as vacancy traps, as in bcc Fe, preliminary calculations 
in Ni show that VC and VN clusters may diffuse cooper- 



atively with an effective migration energy similar to that 
for the isolated vacancy. This also raises the possibility of 
enhanced C and N mobility in irradiated alloys and their 
segregation to defect sinks. A vacancy can bind up to 
two C atoms and up to six N atoms in austenite (or four 
in Ni), although the additional binding energy reduces 
significantly above two. Covalent bonding was observed 
between solutes in a vacancy but did not lead to any en- 
hanced stability, as seen in bcc Fe. Both C and N show 
a strong preference for sites near the surface of vacancy 
clusters and the binding increases with cluster size, sug- 
gesting that they will decorate the surface of voids and 
gas bubbles, when mobile. A binding energy of 0.1 eV 
was observed to a [001] SI dumbbell in austenite and Ni, 
which we would expect to increase with interstitial clus- 
ter size, as in bcc Fe, resulting in Cottrell atmospheres 
of C and N around dislocations and grain boundaries in 
austenitic alloys. 

Along with previous work, these results provide a 
complete database that would allow realistic Fe-Cr-Ni 
austenitic alloy systems to be modelled using higher-level 
techniques, such as molecular dynamics using empirical 
potentials and kinetic Monte Carlo simulations. As such, 
they play a critical role in a multi-scale modelling ap- 
proach to study the microstructural evolution of these 
materials under irradiation in typical nuclear environ- 
ments. 
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